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14.  ABSTRACT 

This  study  is  intended  to  investigate  the  physiological  effects  of  anti-angiogenic  and  anti-tumor  treatment  at  the  microcirculatory  level  with  the  intent  to 
better  understand  how  the  tumor  adapts  to  these  treatments,  with  the  aim  of  formulating  more  effective  treatment  modalities  based  off  this 
knowledge.  We  have  recently  had  the  opportunity  to  form  a  collaboration  with  a  research  group  investigating  the  effects  of  microbeam  radiation 
therapy  (MRT).  Previous  studies  have  shown  that  MRT  induces  a  strong  anti-tumor  effect  while  sparing  normal  tissue.  The  mechanisms  of  this  tissue¬ 
sparing  effect  are  poorly  understood,  however.  Investigation  of  MRT  in  the  murine  window  chamber  model  will  not  only  elucidate  the  mechanisms  of 
MRT,  but  will  also  provide  further  insight  into  natural  tumor  progression  and  how  the  tumor  may  respond  to  vascular  disrupting  agents  in  general.  We 
have  employed  a  number  of  advanced  techniques  in  the  window  chamber  model  to  optically  observe  natural  and  induced  changes  in  tumor  physiology. 
These  advanced  imaging  modalities  have  the  advantage  of  being  able  to  facilitate  longitudinal,  in  vivo  investigation  into  the  parameters  of  interest. 
These  include  Doppler  Optical  Coherence  Tomography  for  the  measurement  of  blood  flow  velocity,  Hyperspectral  Imaging  for  the  measurement  of 
hemoglobin  saturation  and  hematocrit,  and  Confocal  Imaging  for  the  measurement  of  3-D  vascular  architecture  and  perfusion.  Additionally,  novel 
methods  have  been  developed  to  molecularly  probe  the  tumor  environment.  These  include  Oxygen  Sensitive  Nanoparticles  for  the  mapping  of  p02  and 
the  use  of  Engineered  Cell  Lines  which  express  green  fluorescent  protein  under  hypoxia.  Together,  these  techniques  will  allow  us  to  better  understand 
and  model  the  process  of  tumor  growth  and  response  to  MRT,  eventually  leading  to  a  better  understanding  of  how  such  changes  can  be  exploited  to 
increase  the  efficacy  of  combinational  therapies. 
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INTRODUCTION: 


This  study  is  intended  to  investigate  the  physiological  effects  of  anti-angiogenic  and  anti-tumor 
treatment  at  the  microcirculatory  level  with  the  intent  to  better  understand  how  the  tumor  adapts  to 
these  treatments,  with  the  aim  of  formulating  more  effective  treatment  modalities  based  off  this 
knowledge.  We  have  recently  had  the  opportunity  to  form  a  collaboration  with  a  research  group 
investigating  the  effects  of  microbeam  radiation  therapy  (MRT).  Previous  studies  have  shown  that  MRT 
induces  a  strong  anti-tumor  effect  while  sparing  normal  tissue.  The  mechanisms  of  this  tissue-sparing 
effect  are  poorly  understood,  however.  Investigation  of  MRT  in  the  murine  window  chamber  model  will 
not  only  elucidate  the  mechanisms  of  MRT,  but  will  also  provide  further  insight  into  natural  tumor 
progression  and  how  the  tumor  may  respond  to  vascular  disrupting  agents  in  general.  We  have 
employed  a  number  of  advanced  techniques  in  the  window  chamber  model  to  optically  observe  natural 
and  induced  changes  in  tumor  physiology.  These  advanced  imaging  modalities  have  the  advantage  of 
being  able  to  facilitate  longitudinal,  in  vivo  investigation  into  the  parameters  of  interest.  These  include 
Doppler  Optical  Coherence  Tomography  for  the  measurement  of  blood  flow  velocity,  Hyperspectral 
Imaging  for  the  measurement  of  hemoglobin  saturation  and  hematocrit,  and  Confocal  Imaging  for  the 
measurement  of  3-D  vascular  architecture  and  perfusion.  Additionally,  novel  methods  have  been 
developed  to  molecularly  probe  the  tumor  environment.  These  include  Oxygen  Sensitive  Nanoparticles 
for  the  mapping  of  p02  and  the  use  of  Engineered  Cell  Lines  which  express  green  fluorescent  protein 
under  hypoxia.  Together,  these  techniques  will  allow  us  to  better  understand  and  model  the  process  of 
tumor  growth  and  response  to  MRT,  eventually  leading  to  a  better  understanding  of  how  such  changes 
can  be  exploited  to  increase  the  efficacy  of  combinational  therapies. 


BODY: 

The  grant  application's  Statement  of  Work  listed  a  number  of  goals  within  Task  1  to  be  completed  within 
the  first  16  months  of  work.  These  goals  largely  consisted  of  the  development  and  refinement  of 
imaging  modalities  for  the  optical  interrogation  of  tumors  grown  in  a  window  chamber.  While  the 
essence  of  this  project  has  remained  consistent,  we  have  elected  to  employ  the  dorsal  window  chamber 
model  (rather  than  the  mammary  window  chamber  model)  for  the  initial  stages  of  the  project.  While 
the  mammary  window  chamber  model  has  the  advantage  of  allowing  for  orthotopic  tumor 
implantation,  we  felt  that  the  benefits  of  the  dorsal  model  outweighed  this  advantage.  The  ability  to 
use  trans-illumination  to  image  the  dorsal  window  chamber  provides  more  accurate  results  when  using 
the  hyperspectral  imaging  technique,  as  it  does  not  require  complicated  modeling  of  depth-dependant 
tissue  reflectance  and  scattering.  Additionally,  many  of  the  imaging  modalities  would  require  significant 
modification  to  facilitate  the  mammary  model.  For  our  initial  investigation  into  unmodulated  tumor 
progression,  we  felt  that  the  dorsal  model  would  provide  cleaner,  more  easily  processed  data  that 
would  eventually  facilitate  adaptation  to  the  mammary  model,  once  basic  modulations  of  physiological 
parameters  were  thoroughly  modeled  and  understood. 
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Since  submission  of  the  last  annual  report,  we  have  made  significant  refinements  to  the  image 
processing  algorithms  used  to  automatically  quantify  vascular  parameters.  These  algorithms  constitute 
an  invaluable  tool  by  which  we  are  capable  of  quickly  and  objectively  quantifying  induced  changes  in 
vascular  structure  and  functionality.  Variables  that  this  algorithm  is  capable  of  quantifying  include 
vascular  length  density,  total  branchpoints,  tortuosity,  and  vascular  diameter.  Additionally,  we  have 
developed  code  that  is  able  to  process  video-rate  sequences  of  blood  flow  images  and  provide 
quantitative  maps  of  flow  direction  and  velocity.  This  advancement  (in  conjunction  with  other 
techniques  for  measuring  hemoglobin  saturation  and  HIF-1  expression)  is  particularly  useful  in  helping 
us  to  understand  how  blood-flow  patterns  affect  the  local  tumor  microenvironment.  We  expect  to 
publish  data  derived  from  this  algorithm  within  the  next  few  months. 

The  recent  development  of  oxygen  sensitive  nanoparticles  through  collaboration  with  Cassandra  Fraser 
at  the  University  of  Virginia  has  provided  us  with  an  additional  technique  for  optically  measuring 
extravascular  PO21.  We  have  already  employed  this  tool  in  a  preliminary  study  of  oxygen  transport  and 
hemodynamics.  We  will  continue  to  employ  these  nanoparticles  throughout  the  rest  of  this  study. 

Finally,  we  wish  to  report  an  exciting  opportunity  that  had  arisen  soon  after  the  submission  of  the  last 
annual  report.  We  have  been  provided  with  the  opportunity  to  collaborate  with  a  group  at  the 
University  of  North  Carolina  -  Chapel  Hill  that  is  studying  the  effects  of  microbeam  radiation  therapy. 
Microbeam  radiation  therapy  (MRT)  is  a  unique  form  of  radiotherapy  that  has  shown  a  marked  tumor- 
specific  effect.  Compared  to  conventional  radiation  forms,  MRT  produces  microscopic,  spatially-discrete 
radiation  patterns  at  an  ultrahigh  dose  rate.  The  tissue-sparing  property  of  this  unique  treatment 
modality  is  possibly  facilitated  by  efficient  normal  vessel  repair  mechanisms,  contrary  to  the 
catastrophic  disruption  of  poorly  regulated  tumor-associated  vasculature.  However,  the  non-local 
response  characteristics  associated  with  this  treatment  are  not  clearly  understood.  Therefore,  we  have 
begun  work  to  investigate  this  treatment  in  the  window  chamber  model  in  order  to  observe  and 
spatially  quantify  treatment  response  in  terms  of  modulation  of  vascular  structure  and  functionality.  In  a 
preliminary  study,  over  a  week-long  time  course,  we  observed  induced  changes  in  tumor-associated 
vasculature  after  the  MRT  treatment.  Time  course  images  and  analysis  presented  here  suggest  a 
treatment-induced  modulation  of  the  tumor  vasculature  with  regard  to  both  structure  and  function.  We 
will  continue  this  work  throughout  the  next  year,  employing  all  of  the  automated  image  processing 
techniques  that  were  described  previously.  We  expect  that  this  data  will  yield  invaluable  information  on 
how  MRT  is  able  to  specifically  target  the  tumor.  Furthermore,  it  will  provide  additional  information  on 
how  tumor  vasculature  naturally  adapts  its  vascular  delivery  line  to  compensate  for  inefficiencies  in 
function  do  to  poorly  regulated  angiogenesis.  Beyond  investigation  of  the  effects  of  MRT,  these  data 
will  provide  insight  into  the  mechanisms  and  shortcomings  of  traditional  anti-vascular  treatments. 

Below  are  presented  the  tasks  originally  proposed  to  be  completed  in  the  first  16  months  of  the  project, 
along  with  discussions  of  accomplishments  or  alterations  to  these  tasks: 

SOW  -  Task  la:  Learn  mammary  window  chamber  surgical  technique  (month  1-2) 
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As  discussed  above,  we  have  decided  to  employ  dorsal  window  chamber  and  fiber  probe  based 
techniques,  rather  than  the  mammary  window  chamber  procedure.  I  have  since  become  proficient  at 
both  of  these  techniques  and  have  applied  them  in  a  number  of  preliminary  studies. 

SOW  -  Task  lb:  Establish  method  for  extracting  fluorescein  concentrations  in  extravascular 
space  using  published  data  by  Dr.  Gregory  Palmer.  Palmer  will  provide  guidance,  (month  1-2) 

This  method  has  been  established  and  will  be  employed  in  future  drug  kinetics  studies. 

SOW  -  Task  lc:  Hemoglobin  saturation  extraction  software  will  need  to  be  adapted  to  mammary 
window  chamber  measurements,  (month  1-2) 


Figure  1,  showing  the  original  hemoglobin  saturation 
reconstruction:  Hemoglobin  saturation  is  represented  as 
grayscale  values,  independent  of  total  hemoglobin  content 
within  the  vessels.  This  representation  has  the  drawback  of 
not  being  able  to  represent  important  variations  in  individual 
vessels'  hemoglobin  (and  thus  oxygen)  delivering  capacity. 


Figure  2,  showing  the  same  tumor,  but  incorporating 
improvements  made  to  the  hyperspectral  processing 
algorithm:  Total  hemoglobin  is  factored  into  the  new 
algorithm  and  is  represented  as  vascular  brightness. 
Hemoglobin  saturation  is  represented  on  a  blue-red  color 
axis,  with  blue  and  red  representing  fully  deoxygenated  and 
oxygenated  blood,  respectively.  The  tumor  appears  in 
white,  while  hypoxia-indicating  GFP  is  represented  by  its 
green  intensity  within  the  tumor. 


Since  the  dorsal  window  chamber  has  replaced  the  mammary  window  chamber  in  this  aim  of  the  study, 
no  physical  modifications  to  the  hyperspectral  system  were  necessary.  I  have,  however,  made 
improvements  to  the  computational  reconstruction  algorithm  (Figs  1&2),  and  we  are  now  able  to  extract 
total  hemoglobin  values  from  our  data.  Total  hemoglobin  will  be  incorporated  into  all  of  our  window 
chamber  studies  from  this  point  forward. 


SOW -Task  Id:  Practice  imaging  modalities  and  collect  preliminary  data,  (month  3) 


This  task  has  been  completed  and  has  progressed  towards  the  more  detailed  analysis  of  Task  le. 
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SOW  -  Task  le:  Collect  HIF-1, 
hemoglobin  saturation ,  red 
blood  cell  flux,  drug  delivery, 
and  vascular  structure  data 
forSAl.  (month  4-15) 
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Figure  3,  showing  recently  published  data  regarding  longitudinal  variations  in 
important  physiological  parameters  over  a  36  hour  time  course:  Values  at  each  time 
point  were  averaged  across  all  regions  of  interest.  Hemoglobin  saturation,  flow 
velocity,  vessel  diameter,  redox  ratio,  maximum  velocity,  and  shear  rate  were  all 
shown  to  significantly  change  with  time  (p<0.05). 


These  parameters  (minus  the 
HIF-1  reporter),  along  with 
redox  ratio,  were 
investigated  in  two  recent 
papers2'3.  HIF-1  (as  reported 
by  GFP)  presented  an 
obstacle  to  obtaining 
accurate  redox  ratio  values, 
as  the  GFP  spectrum 

overlaps  with  and  overpowers  the  weak  FAD 
fluorescence  signal  (redox  ratio  is  measured  as  the 
ratio  of  FAD  to  NADH).  Since  this  particular 
timecourse  investigated  shortterm  physiological 
changes  (every  6  hours  for  36  hours),  we  felt  that 
redox  ratio  was  a  better  parameter  to  track,  since 
the  long  half-life  of  GFP  might  produce  misleading 
quantities  over  such  a  short  timeframe.  The  data 
that  we  collected  yielded  interesting  correlations 
among  all  of  these  parameters  (Fig  3).  In  a 
separate  paper,  the  oxygen-sensitive  nanoparticles 
were  employed  to  investigate  oxygenation  and 
hemodynamic  properties  in  tumors  (Fig  4).  This 


Figure  4:  The  data  presented  above  are  from  a  recently 
published  study  (Palmer  etal,  I  BO  2010)  in  which 
oxygen  sensitive  nano  particles  probes  were  used  to 
measure  oxygen  tension  in  within  a  tumor  and 
surrounding  tissue.  The  images  show  (a)  the  tumor 
(boxed)  under  bright  field  illumination,  (b)  oxygen 
tension  in  the  percent  of  oxygen,  (c)  hemoglobin  oxygen 
saturation,  (d)  GFP-H1F- 1  fluorescence,  and  (e) 
constituti  vely  expressed  RFP  fluorescence.  The 
hemoglobin  saturation  (c)  is  shown  with  percent  oxygen 
saturation  on  a  color  scale  whose  brightness  is 
modulated  by  the  total  hemoglobin  content(thus,  well- 
vascularized  regions  appearbrighter).  A  combined  plot 
(f)  shows  the  hemoglobin  oxygen  saturation  projected 
onto  a  3-D  surface  whose  vertical  displacement 
corresponds  to  increased  tissue  oxygenation.  A  1mm 
scale  bar  is  shown  in  (a). 


paper  demonstrates  interesting  correlations 
between  p02  and  vascular  function,  and  will 
establish  a  model  by  which  oxygen  delivery  to  the 
tumor  may  be  better  understood. 


SOW  -  Task  If:  Perform  statistical  analysis  on  data 
and  formulate  combinational  treatment  schedules 
for  future  SAs.  These  results  will  be  reviewed  by  a 
statistician,  (month  15-16) 
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The  project  has  progressed 
significantly  ahead  of 
schedule  with  regard  to  this 
particular  task.  We  have 
already  performed  some 
statistical  analysis  on  our 
observations,  and  these 
results  were  published  in  a 
recent  article  (Fig  5)3. 
Furthermore,  we  have 
mostly  finalized  our 
advanced  image  processing 

algorisms  for  the  automated  segmentation  and  analysis  of  vascular  architecture  (Fig  6).  The  extraction 
of  total  hemoglobin  content  was  one  of  the  first  products  of  this  effort,  and  we  are  well  on  our  way 
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Figure  5,  showing  published  statistical  analysis  of  observed  parameters:  Results 
indicated  a  positive  correlation  between  blood  flow,  hemoglobin  saturation,  and 
metabolic  demand  in  tumors.  But  correlations  are  weak,  suggesting  no  single 
biomarker  can  fully  explain  the  behavior  of  another,  and  further  highlighting  the 
need  for  sophisticated  models  of  tumor  growth  and  adaptation. 


Figure  6  shows  a  sample  of  the  results  from  the 
automated  algorithm  used  to  quantify  vascular 
characteristics.  In  the  image  to  the  left,  vessels  are 
identified  and  highlighted  in  red  within  the  region  of 
interested  identified  by  the  green  border.  The  right 
image  shows  the  transformation  of  these  vessel 
outline  data  into  vascular  parameters:  Vessel 
centerlines  are  shown  in  red,  while  blue  points  show 
vessel  branch  points.  The  green  lines  represent 
straight-line  distances  between  branch  points,  from 
which  vessel  tortuosity  can  be  calculated.  Vessel 
diameter  and  vascular  length  density  can  also  be 
calculated  from  these  data. 


towards  developing  an 
all  encompassing 
algorithm  to  quantify 
mean  vessel  density, 
tortuosity,  and  mean 
vascular  diameter  from 
our  images.  Such  an 
algorithm  will  allow  us 
to  automatically  extract 
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unbiased  values  for  these  parameters,  circumventing  the  subjective  and  time-consuming  task  of 
performing  manual  image  analysis. 


Vascular  Parameters 
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Figure  7  shows  results  from  the 
automated  algorithm  used  to 
quantify  vascular 
characteristics.  The  blue  lines 
represent  changes  in  the  low- 
dose  microbeam  animals  over 
time.  Widefield  radiation  and 
mock  treatment  controls  are 
shown  in  red  and  green, 
respectively.  While  these 
preliminary  results  showed 
interesting  qualitative  trends 
and  helped  to  validate  the 
parameter  quantification 
algorithms,  actual  values  for 
these  parameters  yielded 
uncertain  results.  We  expect  to 
see  more  consistent  and 
significant  results  with  a  larger 
number  of  high-dose 
treatments,  however. 
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Figure  8  shows  a  sample  of  the  hemoglobin 
saturation/HIF-1  images  from  the  low-dose 
microbeam  experiment.  Hemoglobin  saturation  is 
representes  as  per  the  colorscale  shown.  HIF-1 
expression  (GFP)  is  represented  in  green.  Low-dose 
treatment  shows  interesting  qualitative  results. 
Clearly  there  is  an  inhibition  in  tumor  growth  and 
vascular  proliferation.  These  characteristics  will  be 
examined  more  thoroughly  with  the  high-dose 
treatments. 


SOW  -  Task  2a:  Collect  HIF-1 ,  hemoglobin 
saturation ,  red  blood  cell  flux,  drug  delivery,  and 
vascular  structure  data  for  SA  2  (17-26) 


As  stated  previously,  focus  on  treatment 
investigation  has  shifted  to  microbeam  radiation 


treatment.  We  have  collected  preliminary  data  investigating  microbeam  patterns  at  lower  doses 


(approximately  10  Gy)  (Fig  7  &  8),  however,  high-dose  (>  25  Gy)  treatments  were  delayed  by  technical 
difficulties  with  the  microbeam  generator.  These  issues  were  resolved  in  late  2010,  and  we  have  begun 


collecting  data  regarding  this  treatment. 
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KEY  RESEARCH  ACCOMPLISHMENTS: 


•  Developed  appropriate  animal  protocols 

•  Improved  hyperspectral  analysis  to  extract  total  hemoglobin  content 

•  Developed  an  algorithm  for  incorporating  multiple  physiological  parameters  into  a  single  image 
for  better  visual  analysis 

•  Developed  methods  for  incorporating  multiple  imaging  modalities  into  a  single  microscope 
system 

•  Collected  data  over  a  36  hour  time  course  which  showed  significant  fluctuations  in  all 
physiological  parameters  over  that  timeframe 

•  Showed  a  statistically  significant  correlation  among  a  number  of  physiological  parameters, 
suggesting  a  multiparametric  model  of  tumor  progression  is  attainable 

•  Investigated  hemodynamic  properties  of  tumors  and  how  these  correspond  to  changes  in 
oxygen  delivery  as  assessed  by  oxygen-sensitive  nanoparticles 

•  Developed  an  algorithm  for  quantification  of  blood  flow  velocity  and  direction  in  video-rate 
image  sequences 

•  Performed  preliminary  investigation  into  low-dose  microbeam  radiation  effects 

•  Began  investigation  of  high-dose  microbeam  radiation  effects 


REPORTABLE  OUTCOMES: 

We  have  published  three  papers  that  demonstrate  significant  progress  towards  our  goal  of  analyzing 
treatment-induced  physiological  changes  and  demonstrating  refined  understanding  of  treatment- 
induced  effects  upon  vascular  structure  and  function. 


Flow-Induced  Pixel  Intensity  Fluctuations 


Right/Left  Movement 


Op/Down  Movei 


nt 


Relative  Blood  Blow  Velocity 


Figure  9  shows  results  from  our  flow 
quantification  algorithm.  A  raw  frame 
from  a  video  image  of  RBC  movement 
is  shown  in  the  top-left.  From  this 
image  sequence,  flow-induced  pixel 
intensity  fluctuations  can  be 
quantified,  and  from  these  data  are 
derived  values  for  left/right 
movement,  up/down  movement,  and 
total  velocity  for  every  pixel  in  the 
video  sequence. 
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CONCLUSION: 


In  the  second  year  of  work  we  have  continued  to  develop  optimized  procedures  for  the  analysis  of  key 
physiological  parameters  using  optical  methods.  Our  results  from  video  analysis  of  red  blood  cell 
movement  are  particularly  interesting  (Fig  9),  and  we  anticipate  incorporating  this  analysis  into  our 
study  and  publishing  this  processing  algorithm  within  the  next  few  months.  In  addition  to  providing  an 
invaluable  tool  for  analysis  of  our  data,  we  expect  that  this  algorithm  will  find  useful  application  in  a 
number  of  instances  of  video  analysis  of  hemodynamics  by  both  our  lab  and  others.  We  have  begun 
investigation  of  the  effects  of  microbeam  radiation  therapy  on  vascular  function,  and  have  already 
completed  a  preliminary  analysis  of  low-dose  treatments.  These  results  provided  and  excellent  model 
upon  which  to  test  our  algorithms.  We  are  now  ready  to  proceed  to  a  large  scale  investigation  of  MRT, 
and  have  recently  begun  a  number  of  animal  experiments. 
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A  dual-emissive-materials  design  concept  enables 
tumour  hypoxia  imaging 

Guoqing  Zhang1,  Gregory  M.  Palmer2,  Mark  W.  Dewhirst2  and  Cassandra  L.  Fraser1* 


Luminescent  materials  are  widely  used  for  imaging  and  sensing 
owing  to  their  high  sensitivity,  rapid  response  and  facile 
detection  by  many  optical  technologies1.  Typically  materials 
must  be  chemically  tailored  to  achieve  intense,  photostable 
fluorescence,  oxygen-sensitive  phosphorescence  or  dual  emis¬ 
sion  for  ratiometric  sensing,  often  by  blending  two  dyes  in 
a  matrix.  Dual-emissive  materials  combining  all  of  these  fea¬ 
tures  in  one  easily  tunable  molecular  platform  are  desirable, 
but  when  fluorescence  and  phosphorescence  originate  from 
the  same  dye,  it  can  be  challenging  to  vary  relative  flu¬ 
orescence/phosphorescence  intensities  for  practical  sensing 
applications.  Heavy-atom  substitution2  alone  increases  phos¬ 
phorescence  by  a  given,  not  variable  amount.  Here,  we  report 
a  strategy  for  modulating  fluorescence/phosphorescence  for 
a  single-component,  dual-emissive,  iodide-substituted  difluo- 
roboron  dibenzoylmethane-polydactic  acid)  (BF2dbm(l)PLA) 
solid-state  sensor  material.  This  is  accomplished  through 
systematic  variation  of  the  PLA  chain  length  in  controlled 
solvent-free  lactide  polymerization3  combined  with  heavy- 
atom  substitution2.  We  demonstrate  the  versatility  of  this 
approach  by  showing  that  films  made  from  low-molecular- 
weight  BF2dbm(l)PLA  with  weak  fluorescence  and  strong  phos¬ 
phorescence  are  promising  as  'turn  onr  sensors  for  aerody¬ 
namics  applications4,  and  that  nanoparticles  fabricated  from 
a  higher-molecular-weight  polymer  with  balanced  fluorescence 
and  phosphorescence  intensities  serve  as  ratiometric  tumour 
hypoxia  imaging  agents. 

Ratiometric  sensing  addresses  challenges  associated  with  sensor 
concentration  and  obviates  the  need  for  specialized  luminescence 
lifetime  instrumentation.  For  example,  Kopelman  and  co-workers 
developed  a  powerful  nanosensor  (PEEBLEs)  for  ratiometric  oxy¬ 
gen  sensing  comprising  a  cyanine  dye  standard  and  an  02 -sensitive 
Pt  porphyrin  phosphor  in  a  sol-gel  matrix  and  demonstrated  its 
utility  for  cell  biology5.  Single-component  dye-polymer  conjugates 
with  both  fluorescence  and  phosphorescence  offer  advantages  over 
three-component  dye/standard/matrix  mixtures.  Dual-emissive 
materials  possess  an  internal  rather  than  external  standard  in 
ratiometric  sensing  schemes,  and  thus,  1:1  fhiorophore/phosphor 
stoichiometry,  resulting  in  greater  sample  homogeneity  and  min¬ 
imal  dye  leaching.  Readily  processable  biomaterials  combined 
with  photostable,  optically  tunable  dyes  that  emit  brightly  even 
in  aqueous  environments  are  beneficial  for  imaging  and  sensing 
in  biomedical  contexts. 

Previously,  we  discovered  that  difluoroboron  dibenzoylmethane- 
poly(lactic  acid)  (BF2dbmPLA;  as  in  Fig.  1,  with  hydrogen  in  place 
of  iodide)  exhibits  unusual,  long-lived  green  room-temperature 
phosphorescence  (RTP)  in  addition  to  intense  blue  fluorescence6. 
Given  that  thermal  decay  pathways  are  often  accessible,  it  is 
rare  to  see  RTP  from  organic  compounds;  typically  organized 


media7,8  or  heavy  atoms9  are  required  for  RTP.  This  feature  adds 
enhanced  capability  to  the  excellent  fluorescence  properties  of 
boron  dyes10,11  and  materials12,13.  Nanoparticles  fabricated  from 
BF2dbmPLA  and  related  analogues  show  great  promise  for  in  vitro 
and  in  vivo  imaging  and  oxygen  sensing14.  Although  BF2dbmPLA 
singlet  and  triplet  emission  are  clearly  distinguished  on  the  basis  of 
wavelength  (AP  =  440  nm,  A.P  =  509  nm)  and  lifetime  (rF  =  2  ns, 
rP  =  170  ms),  the  high  singlet  quantum  yield  (<£F  ~  0.80)  and 
much  weaker  RTP  preclude  the  use  of  this  first-generation  material 
as  a  practical  intensity-based  ratiometric  sensor.  Ideally,  much 
stronger  phosphorescence  intensity  is  required  for  ratiometric 
sensing.  For  insight  into  how  to  achieve  this  goal  and  modulate 
relative  fluorescence  to  phosphorescence  intensities,  we  turn 
to  quantum  theory. 

According  to  perturbation  theory,  for  a  polyatomic  molecule,  the 
perturbation  factor  8  describing  singlet-triplet  mixing  by  first-order 
spin-orbit  coupling  can  be  expressed  as  in  equation  ( 1 )  (ref.  2), 


where  and  are  the  wavefimctions  of  singlet  and  triplet  states, 
respectively,  Hso  is  the  spin-orbit  Hamiltonian  and  Ex  and  E3  are  the 
energy  levels  for  singlet  and  triplet  states.  This  expression  indicates 
that  a  greater  spin-orbit  matrix  value  and  a  smaller  energy  gap 
|£i  —  E3 1  correspond  to  a  larger  8;  that  is,  a  high  probability  of 
electronic  transition  between  singlet  and  triplet  states.  The  spin- 
orbit  matrix  elements  are  highly  sensitive  to  the  atomic  number  of 
the  atoms  in  the  vicinity  of  the  exciton.  Heavy  atoms  either  in  the 
environment  (external),  or  part  of  the  luminescent  molecule  (inter¬ 
nal)  can  enhance  spin-orbit  coupling  and  thus,  increase  the  rate  of 
intersystem  crossing.  Since  the  discovery  of  the  internal  heavy- atom 
effect  by  McClure15  in  1949,  halogens  or  large-atomic-number  met¬ 
als  have  been  incorporated  into  luminescent  molecules  to  enhance 
RTP  intensity  relative  to  fluorescence2.  For  example,  heavy- atom- 
based  RTP  materials  can  improve  the  light- emitting  diode  device16 
efficiency  in  that  triplet  excitons  may  be  the  main  (75%)  excited- 
state  species  on  electron-hole  recombination17.  The  singlet-triplet 
energy  gap,  |Fi  —  E3 1,  on  the  other  hand,  may  be  controlled  by 
boron  dye  loading,  or  in  the  case  of  dye-polymer  conjugates, 
polymer-chain  molecular  weight.  We  have  shown  that  fluorescence 
emission  for  BF2dbmPFA  powders  varies  over  a  64  nm  range  de¬ 
pending  on  BF2dbmPFA  molecular  weight  (3-20  kDa,  507-443  nm 
respectively)3.  High-molecular-weight  samples  exhibit  blueshifted 
emission  similar  to  solution,  whereas  increased  fluorophore- 
fluorophore  interactions  in  low-molecular-weight  materials  corre¬ 
spond  to  lower-energy,  redshifted  fluorescence. 

Phosphorescence,  in  contrast,  showed  little  change  in  energy  or 
relative  intensity  with  BF2dbmPFA  molecular  weight,  suggesting 
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Figure  1 1  Synthesis  and  solid-state  emission  of  BF2dbm(l)PLA  (P1-P3).  a,  Ring-opening  polymerization  of  BF2dbm(l)PLA  (P1-P3).  b-e,  Steady-state 
emission  spectra  of  polymers  P1-P3  as  powders  (b,c)  and  spin-cast  films  (d,e)  under  air  (b,d)  and  under  N2  (c,e).  f,g,  Simple-cast  P1-P3  (right  to  left) 
films  in  vials  under  air  (f)  and  N2  (g)  (ultraviolet  excitation:  A.ex  =  365  nm). 


Table  1 1  Characterization  of  polymers  P1-P3  in  solution. 


Mn 

(Da)* 

PDI ' 

^-abs 

(nm)* 

e 

(M-1  cm"1)* 

A.p 

(nmf 

<2>f1 

tf 

(ns) 

1 

458 

1 

407 

58,000 

441 

0.55 

1.03 

PI 

2,700 

1.11 

406 

33,000 

435 

0.41 

0.95 

P2 

7,300 

1.15 

406 

40,000 

436 

0.44 

0.96 

P3 

17,600 

1.17 

406 

40,700 

435 

0.43 

0.95 

*ln  tetrahydrofuran  versus  polystyrene  standards. 

'  Polydispersity  index  (PDI)  =  Mw/Mn. 

± 

s\  molar  extinction  coefficient.  In  CH2Cl2. 

§ 

Steady-state  fluorescence  spectra  excited  at  369  nm. 

^ Fluorescence  quantum  yields  in  CH2CI2  relative  to  anthracene  in  ethanol. 

II  In  CH2CI2;  excitation  source:  369  nm  light-emitting  diode;  fluorescence  lifetime  fit  to  single¬ 
exponential  decay. 

that  heavy-atom  substitution  may  be  necessary  to  augment  the 
molecular-weight  effect.  On  the  basis  of  these  arguments  and 
previous  observations,  phosphorescence  intensity  is  expected  to 
be  strongest  for  internal  heavy-atom-modified  materials  with 
low  molecular  weights. 

To  test  this  idea,  BF2dbm(I)PLA  samples  with  different  molec¬ 
ular  weights  (PI  =  2,700  Da,  P2  =  7,300  Da  and  P3  =  17,600  Da) 
were  prepared  from  a  BF2dbm(I)OH  initiator,  1,  and  lactide  by 
solvent-free,  tin- catalysed,  controlled  ring-opening  polymerization 
(Fig.  la,  Supplementary  Table  SI).  Polymers  with  low  polydisper¬ 
sity  indices  (PDIs  <  1.2)  were  obtained.  The  optical  properties 
of  P1-P3  were  first  investigated  in  CH2C12  under  air  (Table  1). 
Under  these  conditions,  all  three  polymers  have  nearly  identical 
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Table  2  |  Luminescence  characterization  of  P1-P3  in  the  solid 
state  at  24  °C  under  nitrogen. 


Powder  Film 


*F* 

(nm) 

T 

TF 

(ns) 

^-RTP 

(nm) 

§ 

TRTP 

(ms) 

^F* 

(nm) 

T 

TF 

(ns) 

^-RTP 

(nm) 

§ 

TRTP 

(ms) 

PI 

485 

0.37 

535 

4.06 

458 

0.48 

532 

4.25 

P2 

470 

0.42 

527 

4.39 

445 

0.54 

526 

4.37 

P3 

456 

0.43 

525 

4.50 

438 

0.64 

523 

4.41 

*Steady-state  fluorescence  spectra  excited  at  369  nm. 
i* 

Excitation  source:  369  nm  light-emitting  diode;  fluorescence  lifetime  fit  to  triple¬ 
exponential  decay. 

Jr  Excitation  source:  xenon  flash  lamp  at  405  nm. 

s  Excitation  source:  xenon  flash  lamp  at  405  nm;  RTP  lifetime  fit  to  triple-exponential  decay. 

absorption  (A.abs  =  406  nm)  and  emission  spectra  (A.F  ~  435  nm, 
<2>F  ~  0.4,  rF  ~  0.95  ns)  (see  Supplementary  Figs  SI  and  S2).  As 
expected,  iodide  substitution  results  in  enhanced  intersystem  cross¬ 
ing,  lower  <2>F  and  shorter  rF  compared  with  hydrogen-substituted 
BF2dbmPLA  (ref.  6).  Phosphorescence  is  absent  for  P1-P3  in 
solution;  RTP  is  a  solid-state  effect  for  these  boron  biomaterials. 

For  BF2dbm(I)PLA  powders,  RTP  and  molecular-weight- 
dependent  fluorescence  emission  colour  tuning  are  expected  and 
indeed  observed.  Correspondingly,  two  distinct  emission  bands 
were  observed  in  spectra  for  P2  (470  nm,  527  nm)  and  P3  (456  nm, 
525  nm),  whereas,  for  low- molecular -weight  PI,  a  single  emission 
peak  at  535  nm  was  evident  with  only  a  small  shoulder  at 
~480nm  (Fig.  lc,  Table  2).  The  higher  energy  band  is  attributed 
to  fluorescence  and  the  other  much  stronger,  redshifted  band  is 
RTP  emission.  Under  nitrogen,  the  emission  colours  of  P1-P3 
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Figure  2  |  Oxygen  sensitivity  for  PI  BF2dbm(l)PLA  film,  a,  Emission  spectra  of  the  spin-cast  film  (PI)  under  increasing  oxygen  levels  (indicated  by  the 
arrow,  0-1%)  normalized  to  the  fluorescence  band,  b,  Image  showing  yellow  phosphorescence  emission  under  a  N2  gas  stream  for  a  spin-cast  PI  film 
under  ultraviolet  excitation.  (Yellow  phosphorescence  turns  on  immediately  on  gas  contact.  Blue  green  background:  weak  PI  fluorescence.)  c,  Linear 
relationship  between  oxygen  level  and  the  fluorescence/phosphorescence  intensity  ratio  at  two  fixed  wavelengths  (450  nm  and  525  nm). 


Figure  3  |  Tumour  hypoxia  imaging  with  P2  BF2dbm(l)PLA  nanoparticles,  a-c,  In  vivo  imaging  of  the  breast  cancer  4  T1  mammary  carcinoma  tumour 
region  in  a  mouse  window  chamber  model  showing  the  bright-field  (a)  and  BNP  fluorescence/phosphorescence  ratio  while  breathing  carbogen— 95%  O2 
(b),  room  air— 21%  O2  (c)  and  nitrogen— 0%  O2  (d).  Emission  intensity  was  averaged  from  430  to  480  nm  (fluorescence)  and  530  to  600  nm 
(phosphorescence).  Several  blood  vessels  run  vertically  on  the  left  side  of  the  images  (dark  lines  in  the  bright-field  image;  more  oxygenated  yellow-red 
regions  in  the  fluorescence/phosphorescence  images),  with  the  tumour  comprising  the  region  to  the  right  of  the  vessels  (less-oxygenated  blue  regions  in 
the  fluorescence/phosphorescence  images). 


are  greenish  yellow,  green  and  cyan  respectively  (Fig.  lg).  Just 
as  equation  (1)  predicts,  as  the  energy  gap,  A £,  decreases  (that 
is,  BF2dbm(I)PLA  molecular  weight  decreases),  the  singlet-triplet 
coupling  becomes  stronger  (that  is,  larger  <5),  intersystem  crossing 
is  favoured  and  RTP  intensity  increases. 

Luminescence  lifetime  data  were  also  collected  for  P1-P3  as 
powders  (Table  2).  All  lifetimes  fit  to  triple-exponential  decay 
owing  to  the  heterogeneity  of  the  polymer  matrix  and  possibly 
fluorophore-fluorophore  interactions  too.  The  fluorescence  life¬ 
times  become  shorter  (0.43-0.37  ns)  as  the  polymer  chain  decreases 
(17.6-2.7  kDa)  presumably  owing  to  enhanced  intersystem  crossing 
(Table  1).  Shorter  RTP  and  delayed  fluorescence  lifetimes  (P1-P3: 
4.50-4.06  ms)  with  decreasing  molecular  weight  suggest  smaller 
singlet-triplet  energy  splitting  where  the  thermal  repopulation  from 
triplet  to  singlet  states  is  more  probable3.  Surprisingly,  when  the 
samples  were  exposed  to  air,  long-lived  RTP  was  still  detectable 
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(0.42,  0.55,  0.57  ms  for  P1-P3,  respectively;  see  Supplementary 
Table  S2).  Consistent  with  the  spectral  data  in  Fig.  lb  where  RTP 
bands  remain  visible,  these  results  indicate  a  much  faster  triplet 
decay  rate  compared  with  BF2dbmPLA,  which  is  entirely  quenched 
in  air.  This  suggests  that  radiative  decay  occurs  on  a  comparable 
timescale  to  dynamic  quenching,  which  is  probably  influenced  by 
02  diffusion  rates  in  the  material18. 

Although  the  model  predicts  the  observed  fluorescence/phos- 
phorescence  control  through  heavy-atom  assistance  and  polymer 
molecular  weight,  the  experimental  results  exceeded  our  expec¬ 
tations.  We  discovered  a  tunability  range  that  is  nearly  ideal,  in 
terms  of  control  and  flexibility  in  oxygen  sensor  design.  All  P1-P3 
materials,  regardless  of  molecular  weight  can  be  useful  in  lifetime 
or  time-gated  intensity-based  sensing  platforms.  Those  with  com¬ 
parable  and  readily  detectable  fluorescence  and  phosphorescence 
intensities  (for  example,  P2)  are  beneficial  for  ratiometric  hypoxia 
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imaging,  whereas,  low-molecular-weight  materials  with  weak  to 
negligible  fluorescence  and  strong  phosphorescence  (PI)  serve 
as  ‘turn  on’  sensors  that  light  up  in  anaerobic  or  low- oxygen 
environments.  A  simple  demonstration  of  these  effects  is  presented 
in  Fig.  1,  which  shows  changes  in  colour  and  brightness  for  P1-P3 
films  under  air  (Fig.  If)  versus  nitrogen  (Fig.  lg).  Thus,  with  a 
very  simple  molecular- weight  adjustment,  these  boron  biomateri¬ 
als  are  easily  tailored  for  different  sensor  platforms  and  imaging 
and  detection  schemes. 

To  explore  oxygen  sensitivity  and  processing  effects  in  more 
detail  and  to  demonstrate  the  relevance  and  versatility  of  this 
materials  design  concept  for  practical  applications,  boron  bioma¬ 
terials  were  fabricated  as  thin  films  and  nanoparticles.  Polymers 
P1-P3  were  cast  from  CH2C12  solutions  (5  mg  ml  l,  or  0.3%  w/w) 
onto  glass  substrates,  followed  by  annealing  at  ~  110  °C  for  10  s  to 
improve  film  clarity.  Under  nitrogen,  all  three  films  showed  visibly 
blueshifted  colours  compared  with  the  powders.  The  steady-state 
emission  spectra  under  air  and  nitrogen  are  given  in  Fig.  Id  and  e 
respectively,  where  it  is  noted  that  the  fluorescence  contributions 
are  clearly  larger  than  for  the  bulk  polymer  samples.  Film  fluo¬ 
rescence  lifetimes  are  longer  than  their  powder  counterparts,  sug¬ 
gesting  decreased  intersystem  crossing  perhaps  due  to  diminished 
fluorophore-fluorophore  interactions. 

As  a  starting  point  for  ratiometric  oxygen  sensing,  the  steady- 
state  emission  spectra  of  the  PI  film  were  recorded  under  different 
02  concentrations  (Fig.  2a).  The  ratio  of  the  invariant  fluorescence 
(AF  =  458  nm)  to  the  oxygen-dependent  RTP  (AP  =  532  nm)  steadily 
increases  with  increasing  oxygen  levels  (Fig.  2c).  Up  to  1%  02, 
a  linear  dependency  is  observed  ( R 2  =  0.996).  Beyond  this  point, 
the  fluorescence/phosphorescence  plot  continues  to  rise  up  to 
ambient  levels  (that  is,  21%  02)  but  with  curvature  (Supplementary 
Fig.  S3),  possibly  due  to  complicated  dynamics  in  the  heterogeneous 
polymer  environment  as  previously  reported6.  Figure  2b  shows  the 
oxygen  distribution  on  the  surface  of  a  PI  film  under  ultraviolet 
excitation  (Aex  =  365  nm),  where  the  bright  yellow  streaks  provide 
clear  visualization  of  nitrogen  gas  flow  dynamics. 

The  polymer  P2,  with  balanced  fluorescence  and  phos¬ 
phorescence  emission,  was  also  fabricated  as  boron  nanoparticles 
(BNPs;  98  nm  diameter)  by  nanoprecipitation14.  The  emission 
spectra  (AP  =  450  nm,  ARTP  =  528  nm;  Supplementary  Fig.  S4), 
lifetimes  (rF  ~  0.45  ns,  rRTP  ~  4.82  ms)  and  oxygen  calibration 
(Supplementary  Fig.  S5)  are  comparable  to  P2  films.  Fortuitously, 
the  optimal  linear  range  for  this  nanosensor  (~0-3%)  is  in  excellent 
correspondence  with  hypoxia  in  biological  contexts.  For  example, 
tumour  hypoxia  ( p02  <  1%;  ref.  19)  is  associated  with  increased 
invasiveness  and  resistance  to  radiation  and  chemotherapy.  Despite 
its  importance  for  cancer  treatment  and  tumour  biology20,  hypoxia 
is  difficult  to  image  with  good  spatial  and  temporal  resolution,  par¬ 
ticularly  in  combination21.  Given  their  small  size,  biocompatibility, 
photostability,  dual- emissive  features  and  high  oxygen  sensitivity, 
BNPs  have  the  potential  to  address  some  of  these  challenges.  To 
test  their  ability  as  ratiometric  tumour  hypoxia  imaging  agents, 
P2  nanoparticles  were  used  in  combination  with  a  mouse  dorsal 
window  chamber  breast  cancer  4T1  mammary  carcinoma  model 
for  hyperspectral  imaging22.  Tissue  oxygen  maps  of  the  window 
region  presented  as  fluorescence/phosphorescence  ratios  during 
carbogen,  air  and  brief  nitrogen  breathing  (95,  21  and  0%  02 
respectively),  show  excellent  contrast  between  the  micro  vascula¬ 
ture  (red)  and  the  tumour  tissue  (blue),  which  remained  hypoxic 
regardless  of  the  breathing  gas  (Fig.  3).  The  signal  was  relatively 
stable  ~l-2  min  after  changing  the  breathing  gas,  which  may  be  a 
function  of  the  biology  of  tissue  oxygen  delivery  and  consumption 
more  so  than  the  inherent  response  of  the  nanosensors.  The  BNP 
tumour  oxygenation  maps  are  complementary  to  existing  optical 
methods  such  as  haemoglobin  saturation  imaging,  which  provides 
vascular  oxygenation22. 


We  have  devised  a  simple  yet  powerful  strategy  for  systematically 
tuning  relative  fluorescence  and  RTP  intensities  in  a  dual- 
emissive  boron  biomaterial.  The  experimental  results  support  the 
proposed  design  concept  involving  manipulation  of  both  spin- 
orbit  coupling  and  singlet-triplet  energy  splitting.  The  wide-range 
fluorescence/phosphorescence  tunability  of  boron  biomaterials 
enables  application  in  lifetime,  time-gated  intensity,  ratiometric 
or  ‘turn  on’  sensing  modes,  and  shows  great  promise  for 
optical  imaging,  oxygen  sensing  and  aerodynamics  applications4. 
A  new  method  for  quantifying  tumour  hypoxia  has  also  been 
demonstrated,  and  can  further  understanding  of  the  relationship 
between  hypoxia,  tumour  progression,  metastasis  and  treatment 
resistance.  This  first  demonstration  of  hypoxia  imaging  bodes  well 
for  using  boron  biomaterials  for  oxygen  sensing  in  other  biomedical 
and  biotechnological  contexts. 

Methods 

Molecular  weights  were  determined  by  gel  permeation  chromatography 
(tetrahydrofuran,  20  °C,  1.0  ml  min-1)  versus  polystyrene  standards  on  a 
Hewlett-Packard  instrument  (series  1100  HPLC)  equipped  with  Polymer 
Laboratories  5  pm  mixed- C  columns  and  connected  to  ultraviolet-visible  and 
refractive  index  (Viscotek  LR  40)  detectors.  A  0.58  correction  factor  was  applied23,24. 
Data  were  processed  with  the  OmniSEC  software  (version  4.2,  Viscotek  Corp). 
Ultraviolet-visible  spectra  were  recorded  on  a  Hewlett-Packard  8453  diode  array 
spectrophotometer.  Photographs  were  taken  in  the  dark  using  a  Canon  PowerShot 
SD600  Digital  Elph  camera  with  the  automatic  setting  (no  flash). 

Powders  were  analysed  as  precipitated.  A  Laurell  Technologies  WS-650S 
spin-coater  was  used  to  cast  polymer  films  for  luminescence  measurements  with  the 
default  set-up  (30  s  at  a  constant  speed  of  4,000  r.p.m.).  Boron  polymer  films  were 
spin-cast  from  CH2C12  solutions  (~2.5%  w/w)  onto  Fischer  Scientific  glass  cover 
slides  (22  x  22  mm)  at  a  spin  speed  of  4,000  r.p.m.  The  coated  slides  were  cut  into 
~3  mm  x  22  mm  strips,  placed  in  a  transparent  glass  vial  under  a  N2  atmosphere 
and  were  sealed  with  Teflon  caps  for  measurements.  (Note:  both  borosilicate  glass 
cover  slides  and  glass  vials  are  optically  inactive  using  an  excitation  wavelength 
of  greater  than  368  nm.) 

Steady-state  fluorescence  emission  spectra  were  recorded  on  a  Horiba 
Fluorolog-3  Model  FL3-22  spectrofluorometer  (double-grating  excitation  and 
double-grating  emission  monochromators).  For  all  luminescence  measurements, 
a  Tiffen  55  mm  0.9  neutral  density  filter  was  placed  in  the  path  of  the  excitation 
source  to  protect  the  sample  from  photobleaching.  RTP  spectra  were  recorded 
with  the  same  instrument  except  that  a  pulsed  xenon  lamp  (A.ex  =  369  nm; 
duration  <  1  ms)  was  used  and  spectra  were  collected  with  a  1  ms  delay  after 
excitation.  Time-correlated  single-photon  counting  fluorescence  lifetime 
measurements  were  carried  out  with  a  NanoLED-370  (369  nm)  excitation  source 
and  a  DataStation  Hub  as  the  single-photon  counting  controller.  Phosphorescence 
lifetimes  were  measured  with  a  500  ns  multi  channel  scalar  card  excited  with  a 
pulsed  xenon  lamp  (A.ex  =  369  nm;  duration  <  1  ms).  Lifetime  data  were  analysed 
with  DataStation  v2.4  software  from  Horiba  Job  in  Yvon. 

Fluorescence  quantum  yields,  <£F,  for  BF2dbm(I)OH  (1)  and  BF2dbm(I)PLA 
(P1-P3)  in  CH2C12  were  calculated  versus  anthracene  in  ethanol  as  a  standard,  as 
previously  described25  using  the  following  values:  ^F  anthracene  =  0.27  (ref.  26), 
rip  ethanol  =  1.360,  CH2C12  =  1.424;  n ^  is  the  refractive  index  (20  °C,  yellow 

sodium  D  line,  589  nm).  Optically  dilute  CH2C12  solutions  of  BF2dbm(I)OH  and 
BF2dbm(I)PLA,  and  ethanol  solutions  of  the  anthracene  standard  were  prepared 
in  1-cm-path-length  quartz  cuvettes,  and  absorbances  (A  <  0.1)  were  recorded  and 
steady-state  emission  spectra  were  obtained  (A.ex  =  350  nm;  emission  integration 
range:  365-700  nm). 

For  oxygen  sensitivity  and  intensity-based  Stern-Volmer  plot  measurements, 
gas  mixtures  of  02  and  N2  of  various  concentrations  were  prepared  by  using 
two  Cole-Parmer  65  mm  flow  meters.  For  each  concentration,  the  pre-mixed  gas 
was  passed  through  a  1-m-path  quartz  fluorometer  cell  equipped  with  a  septum 
screw  top  containing  the  spin-cast  PI  film  for  15  min  and  then  the  spectrum  was 
recorded.  The  oxygen  sensitivity  for  the  P2  nanoparticle  aqueous  suspension  was 
calibrated  similarly  except  that  the  pre-mixed  gas  was  first  evenly  dispensed  into  a 
distilled-water  chamber  before  passing  through  a  quartz  cell  containing  1  ml  of  the 
optically  dilute  (abs  <  0.1)  sample. 

A  window  chamber  was  implanted  using  previously  described  techniques22. 
Briefly,  the  mouse  was  anaesthetized  using  an  intraperitoneal  injection  of 
100  mg  kg-1  ketamine,  10  mg  kg-1  xylazine.  The  skin  on  the  back  of  the  animal 
was  stretched  into  a  metal  frame  to  enable  insertion  of  a  titanium  window  frame, 
having  a  12-mm-diameter  window.  The  front  face  of  the  skin  fold  was  excised,  and 
approximately  20,000  4T1  murine  mammary  carcinoma  cells  were  injected  into  the 
underlying  fascia.  A  cover  glass  was  placed  over  the  open  face  of  the  window,  and 
the  tumour  was  allowed  to  incubate  for  10  days  before  imaging. 

An  approximately  100  pi  suspension  of  P2  nanoparticles  (~1  mg  ml-1) 
was  injected  into  the  space  between  the  window  chamber  and  the  cover  glass. 
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Approximately  5  min  passed  between  the  injection  and  imaging,  which  allowed 
for  stabilization  of  the  signal.  The  animal  was  placed  on  an  upright  fluorescence 
microscope  equipped  with  a  DAPI  excitation  filter,  and  a  liquid-crystal  tunable 
emission  filter  (VariSpec,  Cambridge  Research  and  Instrumentation).  This  enabled 
acquisition  of  an  emission  spectrum,  which  was  acquired  in  10  nm  intervals  from 
430  to  600  nm.  Gas  was  administered  by  a  nose  cone  at  a  rate  of  5 1  min-1,  except 
for  room  air,  during  which  the  gas  was  shut  off.  The  carbogen  breathing  was  begun 
at  the  beginning  of  the  experiment,  with  approximately  5  min  passing  before  the 
initial  acquisition.  The  gas  was  then  switched  to  room  air  for  a  period  of  3  min,  and 
finally  nitrogen  for  30  s.  The  images  shown  in  Fig.  3  were  taken  at  the  end  of  each 
respective  period  of  gas  breathing.  Data  were  processed  by  taking  the  ratio  of  the 
fluorescence  to  phosphorescence  signals,  which  were  defined  as  the  average  signal 
acquired  from  430  to  480  nm  and  530  to  600  nm,  respectively. 
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We  have  combined  hyperspectrai  imaging  with  spectral  domain  optical  coherence  tomography  (SDOCT)  to 
noninvasively  image  changes  in  hemoglobin  saturation,  blood  flow,  microvessel  morphology,  and  sheer  rate 
on  the  vessel  wall  with  tumor  growth.  Changes  in  these  hemodynamic  variables  were  measured  over  24  h  in 
dorsal  skin  fold  window  chamber  tumors.  There  was  a  strong  correlation  between  volumetric  flow  and  he¬ 
moglobin  saturation  (p=0.89,  p  =  9  X  10-6,  N=  15)  and  a  moderate  correlation  between  shear  rate  on  the  ves¬ 
sel  wall  and  hemoglobin  saturation  (p=0.56,  p  =  0.03,  N=  15).  ©  2009  Optical  Society  of  America 
OCIS  codes:  170.4500,  110.4234,  170.2655,  170.3880. 


Hypoxia  is  important  in  regulating  tumor  aggressive¬ 
ness  and  treatment  resistance  and  thus  serves  as  a 
valuable  biomarker  for  tumor  treatment  outcome.  A 
better  understanding  of  tumor  microcirculation 
would  provide  a  basis  for  improved  prognostic  and 
treatment  approaches  for  solid  tumors.  In  tumors, 
formation  of  new  vessels  and  remodeling  and  dropout 
of  existing  vessels  occurs  continuously,  resulting  in 
temporal  variations  in  flow  and  oxygenation  [1].  Dy¬ 
namic  imaging  of  vessel  morphology,  hemoglobin  oxy¬ 
gen  saturation,  blood  flow,  and  shear  rate  on  the  ves¬ 
sel  wall  would  provide  valuable  insight  into  the 
mechanisms  and  distribution  of  cycling  hypoxia  in 
growing  tumors.  For  example,  these  data  will  allow 
for  direct  observation  of  tumor  microvessel  response 
to  changing  hemodynamic  and  metabolic  conditions 
induced  by  cancer  treatments. 

Tumor  hemodynamics  on  the  microvascular  level 
have  traditionally  been  studied  with  invasive  meth¬ 
ods  or  highly  localized  (nonimaging)  methods.  Confo- 
cal  or  multiphoton  microscopy  coupled  with  a  vascu¬ 
lar  contrast  agent  has  been  widely  used  in  animal 
models  to  track  changes  in  vessel  morphology  [2]. 
Blood  flow  in  small  vessels  has  been  limited  to  point 
measurements  using  video-rate  microscopy  and 
single  cell  counting  [1]  or  laser  Doppler  flowmetry 
[3] .  Tumor  oxygenation  has  been  quantified  with  fluo¬ 
rescence  lifetime  imaging  of  a  p02  calibrated  dye  [4] 
or  with  microelectrode  measurements  [3] .  Point  mea¬ 
surements  of  shear  rates  on  vessel  walls  have  also 
been  estimated  with  video  microscopy  [5].  Noninva¬ 
sive  imaging  methods  are  currently  emerging  to  re¬ 
place  these  traditional  techniques.  Speckle  variance 
optical  coherence  tomography  (OCT)  has  recently 
shown  promise  for  monitoring  changes  in  vessel  mor¬ 
phology  with  photodynamic  therapy  in  the  dorsal 
skin  fold  window  chamber  [2] .  Doppler  OCT  imaging 
of  blood  flow  direction  and  velocity  profiles  within 
vessels  [6],  and  hyperspectrai  imaging  of  tumor  he¬ 
moglobin  oxygenation  saturation  [7],  have  also  been 
shown  separately  to  provide  valuable  information  on 
tumor  vessel  function.  However,  a  complete  picture  of 


tumor  hemodynamics  requires  multidimensional  im¬ 
age  data  including  vessel  morphology,  blood  velocity 
profiles  and  direction  of  flow,  and  hemoglobin  oxygen 
saturation. 

We  have  combined  hyperspectrai  imaging  with 
OCT  to  noninvasively  image  changes  in  vascular 
morphology,  hemoglobin  oxygen  saturation,  blood  ve¬ 
locity  and  direction  of  flow,  and  shear  rate  on  the  ves¬ 
sel  wall  with  tumor  growth.  This  combined  micro¬ 
scope  provides  a  wealth  of  information  on  the 
dynamics  of  structural  and  functional  changes  in  tu¬ 
mor  vasculature.  Hyperspectrai  (hemoglobin  oxygen 
saturation)  and  spectral  domain  OCT  (SDOCT) 
(blood  flow  and  vascular  morphology)  arms  were  con¬ 
nected  through  two  separate  baseports  in  an  inverted 
microscope  (Carl  Zeiss  Axiovert  200)  (Fig.  1).  Two- 
dimensional  hyperspectrai  images  were  collected 
with  a  100  W  halogen  lamp  for  transillumination, 
and  detection  was  achieved  with  a  liquid-crystal  tun¬ 
able  filter  (LCTF)  (CRI,  Inc.)  placed  in  front  of  a  DVC 
1412  CCD  camera  (DVC  Company)  [7].  Custom  soft¬ 
ware  was  used  to  tune  the  filter  and  acquire  images 
at  10  nm  increments  between  500  and  610  nm.  Mea¬ 
surements  of  the  dark  offset  and  transmission 
through  a  neutral  density  filter  at  each  wavelength 
were  made  before  each  imaging  session.  Hemoglobin 
saturation  images  were  calculated  by  applying  an  ex¬ 
tension  of  the  Beer-Lambert  law  to  the  wavelength- 
dependent  absorption  at  each  pixel  (assuming  oxy¬ 
genated  and  deoxygenated  hemoglobin  as  the  only 
absorbers),  and  then  solving  for  hemoglobin  satura¬ 
tion  with  linear  least-squares  regression  [7] .  The  sys- 
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Fig.  1.  Schematic  of  the  combined  hyperspectral-SDOCT 
microscope.  L1-L4,  lenses  1-4;  TL,  tube  lens;  BS,  beam 
splitter;  LCTF,  liquid-crystal  tunable  filter. 
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tem,  software,  and  analysis  techniques  have  previ¬ 
ously  been  validated  on  liquid  phantoms  with  an 
accuracy  of  approximately  1%  and  in  vivo  [7] . 

The  common-path  spectral  domain  OCT  arm  con¬ 
sists  of  a  Tiisapphire  laser  source  (Femtolasers)  cen¬ 
tered  at  790  nm  with  a  90  nm  FWHM.  The  interfero- 
gram  is  detected  using  a  custom-made  spectrometer 
with  a  line-scan  CCD  camera  (Atmel,  Aviiva).  The 
OCT  system  is  driven  by  software  that  controls  the 
lateral  scanner  and  performs  data  acquisition  and  ar¬ 
chiving  (Bioptigen,  Inc.).  Three-dimensional  blood 
flow  images  were  collected  with  Doppler  OCT,  which 
measures  phase  changes  due  to  flowing  erythrocytes. 
Phase  changes  are  calculated  from  multiple  A  scans 
collected  at  the  same  position  in  the  volume,  and  can 
be  related  to  flow  velocity  profiles  by  measuring  the 
angle  of  incidence  in  the  three-dimensional  volume. 
Three-dimensional  morphology  was  collected  with 
speckle  variance  OCT,  which  detects  blood  vessels 
from  the  variance  in  the  speckle  pattern  between  se¬ 
quential  B-scan  magnitude  images  [2].  Speckle  vari¬ 
ance  is  advantageous  for  morphology  mapping,  be¬ 
cause  it  is  independent  of  the  angle  between  the 
blood  flow  and  the  incident  beam.  However,  it  is  cur¬ 
rently  not  capable  of  accurately  quantifying  blood  ve¬ 
locity  or  flow  direction. 

Hyperspectral,  Doppler  OCT,  and  speckle  variance 
OCT  images  were  collected  from  4T1  tumors  im¬ 
planted  in  dorsal  skin  fold  window  chambers  in  nude 
mice.  Surgery  and  imaging  were  carried  out  under 
ketamine-xylazine  and  isoflurane  anesthesia,  re¬ 
spectively,  with  the  mice  maintained  at  body  tem¬ 
perature  (conducted  with  institutional  approval  at 
Duke  University).  The  transmission  image  of  the 
window  chamber  morphology  (Fig.  2A)  and  hyper¬ 
spectral  image  of  the  percent  hemoglobin  oxygen¬ 
ation  saturation  in  the  vessels  (Fig.  2C)  were  taken 
with  a  2.5X  objective  (NA=0.12).  Speckle  variance 
OCT  (Fig.  2B)  provides  three-dimensional  vessel 
morphology,  and  Doppler  OCT  (Fig.  2D)  provides 
three-dimensional  vessel  flow  velocities  and  flow  di¬ 
rection.  Speckle  variance  and  Doppler  OCT  images 
were  collected  over  a  1  mm  X  1  mm  area  with  250 
X  125  pixels,  and  1024  pixels  in  the  depth  dimension. 
OCT  images  were  collected  with  a  4x  objective  (NA 
=  0.1)  with  a  1  ms  integration  time  for  each  A  line. 
Eight  repeated  B  scans  were  collected  for  speckle 
variance  OCT,  registered  using  “stackreg”  for  Im- 
ageJ,  and  10  repeated  A  scans  were  colleted  for  Dop¬ 
pler  OCT.  The  flow  profiles  in  Figs.  2E  and  2F  were 
fit  to  a  second-order  polynomial  and  corrected  for  the 
angle  of  incidence  to  provide  velocity  in  millimeters 
per  second.  Maximum  blood  velocity  is  determined 
from  the  peak  of  this  fit,  vessel  diameter  is  deter¬ 
mined  from  the  zero  crossings  of  the  fit,  and  the 
shear  rate  on  the  vessel  wall  is  calculated  from  the 
derivative  of  the  fit  {dvjd r),  assuming  a  Newtonian 
fluid  [8] .  Coregistration  of  speckle  variance  OCT,  hy¬ 
perspectral,  and  Doppler  OCT  images  allow  for  vessel 
morphology,  percent  hemoglobin  oxygen  saturation, 
blood  velocity  and  direction  of  flow,  and  shear  rate  on 
the  vessel  to  be  determined  at  any  point  within  the 
window  chamber. 


Fig.  2.  (Color  online)  Multidimensional  functional  imag¬ 
ing  of  4T1  tumor  vasculature  in  the  dorsal  skin  fold  mouse 
window  chamber.  Transmission  image  of  the  window  cham¬ 
ber  taken  with  the  liquid  crystal  tunable  filter  set  at 
500  nm,  with  a  portion  of  the  tumor  visible  (arrow,  A). 
Scale  bar  is  200  jam.  B,  Speckle  variance  OCT  image  of  ves¬ 
sel  morphology  shown  as  an  en  face  average  intensity  pro¬ 
jection  over  1  mm  depth.  C,  Hyperspectral  image  of  percent 
hemoglobin  oxygen  saturation,  thresholded  for  the 
R-squared  value  of  the  linear  fit  and  total  absorption  value 
at  each  pixel  (background  pixels  are  black).  D,  Doppler 
OCT  image  of  vessel  blood  flow  direction  shown  as  an  en 
face  maximum  intensity  projection  over  depth,  with  red 
vessels  (1)  flowing  toward  the  top  of  the  image  and  blue 
vessels  (2)  flowing  toward  the  bottom  of  the  image.  The 
flow  profiles  at  points  1  (E)  and  2  (F)  were  fit  to  a  second- 
order  polynomial  and  corrected  for  the  angle  of  incidence  to 
provide  velocity  in  millimeters  per  second. 

Changes  in  hemoglobin  saturation,  vessel  morphol¬ 
ogy,  and  blood  flow  with  tumor  growth  were  mea¬ 
sured  in  a  14-day-old  4T1  tumor  in  the  dorsal  skin 
fold  window  chamber  (Fig.  3).  Hyperspectral  and 
Doppler  OCT  data  sets  were  collected  every  6  h  for 
24  h.  The  blood  velocity  and  hemoglobin  saturation 
as  a  function  of  time  are  shown  for  three  vessel  cross 
sections  of  interest  (VOIs)  in  Fig.  3.  The  measured 
hemoglobin  saturation  values  fall  within  the  ex¬ 
pected  range  [9,10].  Note  that  the  p02  of  peripheral 
arteries  in  both  tumors  [9]  and  normal  tissues  [10] 
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Fig.  3.  (Color  online)  Time  series  data  taken  from  a  4T1 
window  chamber  tumor  using  the  combined  hyperspectral- 
OCT  microscope.  The  transmission  image  (A)  is  3  mm 
(height)  X  2  mm  (width)  and  delineates  the  three  VOIs  for 
which  maximum  blood  velocity  (B),  diameter  (C),  and  he¬ 
moglobin  saturation  (D)  are  plotted  (open  circles,  VOI  1; 
open  triangles,  VOI  2;  open  squares,  VOI  3). 
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can  be  substantially  less  than  95%  saturated  owing 
to  loss  of  oxygen  from  the  peripheral  arterial  tree 
(longitudinal  oxygen  gradient).  The  extracted  maxi¬ 
mum  velocities  and  vessel  diameters  are  in  agree¬ 
ment  with  previously  published  values  from  the  win¬ 
dow  chamber  model  [11].  Variations  in  red  blood  cell 
flux  and  p 02  [1],  hemoglobin  oxygen  saturation  [7], 
and  diameter  [9]  have  previously  been  observed  in 
tumors  over  time  scales  ranging  from  one  hour  to 
several  days. 

The  relationship  between  the  maximum  velocity, 
diameter,  flow  rate,  shear  rate,  and  hemoglobin  oxy¬ 
gen  saturation  for  the  same  three  vessels  is  shown  in 
Fig.  4.  Spearman  rank  correlation  analysis  of  the 
three  vessels  over  five  time  points  revealed  a  strong 
correlation  (p=0.94,  p  =  0.018,  N=  5)  between  vessel 
hemoglobin  saturation  and  maximum  vessel  velocity 
for  VOI  2  (Fig.  4B)  and  no  significant  correlation  for 
VOI 1  or  VOI  3  (p  >  0.05,  Figs.  4A  and  4C).  There  was 
no  correlation  between  vessel  diameter  and  maxi¬ 
mum  velocity  or  between  vessel  diameter  and  hemo¬ 
globin  saturation  for  any  of  the  VOIs  (p>0.05).  The 
flow  rate  [Q  =  7t*(Z)/2)2*  Vmax/2(mm3/s)]  was  calcu¬ 
lated  from  the  vessel  diameter  (. D )  and  maximum 
vessel  velocity  (Vmax).  For  all  three  VOIs  and  five 
time  points  grouped  together,  there  was  a  strong  cor¬ 
relation  between  velocity  and  hemoglobin  saturation 
(Fig.  4D,  p=  0.78,  p  =  0.0006,  iV=  15),  and  between  flow 
rate  and  hemoglobin  saturation  (Fig.  4E,  p=  0.89,  p 
=  9XlO-6,  iV=15),  and  a  moderate  correlation  be¬ 
tween  shear  rate  on  the  vessel  wall  and  hemoglobin 
saturation  (Fig.  4F,  p= 0.56,  p  =  0.03,  iV=15).  Shear 
rate  values  are  in  qualitative  agreement  with  previ¬ 
ous  studies  employing  different  techniques  [5] . 

Previous  studies  [1,4]  indicate  that  vessel  remodel¬ 
ing  with  tumor  growth  alters  flow  distributions,  and 


Fig.  4.  (Color  online)  Scatter  plots  of  maximum  velocity 
vs.  hemoglobin  saturation  (O)  and  diameter  vs.  hemoglobin 
saturation  (A)  for  VOI  1  (A),  VOI  2  (B),  and  VOI  3  (C). 
There  was  a  strong  correlation  (B)  between  the  maximum 
velocity  and  hemoglobin  saturation  for  VOI  2  (p= 0.94,  p 
=  0.018,  N=  5),  and  no  correlation  for  any  other  VOI  (p 
>0.05).  However,  when  all  VOI  (O,  VOI  1;  A,  VOI  2;  □, 
VOI  3)  and  time  points  were  grouped  together  (D),  there 
was  a  strong  correlation  between  maximum  velocity  and 
hemoglobin  saturation  (p= 0.78,  p  =  0.0006,  N=15).  There 
was  also  a  strong  correlation  (E)  between  flow  (vessel  cross 
sectional  area  times  \  maximum  velocity)  and  hemoglobin 
saturation  (p=0.89,p  =  9x  10-6,  iV=15)  and  a  moderate  cor¬ 
relation  (F)  between  shear  rate  on  the  vessel  wall  and  he¬ 
moglobin  saturation  (p = 0.56,  p-  0.03,  N- 15)  when  all  VOI 
and  time  points  were  grouped  together. 


the  relationship  between  hemoglobin  oxygen  satura¬ 
tion  and  blood  flow  in  tumors  is  heterogeneous.  How¬ 
ever,  the  techniques  employed  in  previous  studies 
have  a  significant  effect  on  the  results.  The  dual  slit 
technique  suffers  from  inaccurate  estimations  of  the 
length  traveled  by  red  blood  cells,  because  the  length 
of  a  three-dimensional  vessel  is  projected  onto  a  two- 
dimensional  image,  resulting  in  imprecise  velocity 
values.  Red  blood  cell  flux  measurements  are  limited 
to  vessels  smaller  than  —30  pm  in  diameter,  because 
individual  red  blood  cells  are  difficult  to  visualize  in 
larger  vessels.  Also,  previous  studies  measured  intra¬ 
vascular  and  extravascular  p 02,  rather  than  hemo¬ 
globin  saturation.  The  techniques  demonstrated  in 
the  current  study  are  advantageous,  because  flow  ve¬ 
locity  profiles  can  be  accurately  measured  from  a 
wide  range  of  vessel  diameters,  and  morphology, 
blood  oxygen  saturation,  and  blood  velocity  can  be  ac¬ 
curately  determined  without  the  use  of  contrast 
agents. 

We  have  combined  hyperspectral  microscopy  with 
OCT  to  create  a  tool  that  can  dynamically  and  non- 
invasively  monitor  vessel  structure  and  function. 
These  combined  measurements  will  allow  for  unprec¬ 
edented  insight  into  the  relationship  among  vessel 
structure,  hemoglobin  saturation,  blood  flow,  and  hy¬ 
poxic  episodes  in  tumors.  This  research  instrument 
could  ultimately  lead  to  a  more  detailed  understand¬ 
ing  of  the  structural  and  functional  properties  of  tu¬ 
mor  vasculature  and  potentially  become  a  valuable 
tool  for  screening  tumor  therapies  in  animal  models. 
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Abstract.  An  important  feature  of  tumor  hypoxia  is  its  temporal  insta¬ 
bility,  or  "cycling  hypoxia."  The  primary  consequence  of  cycling  hy¬ 
poxia  is  increased  tumor  aggressiveness  and  treatment  resistance  be¬ 
yond  that  of  chronic  hypoxia.  Longitudinal  imaging  of  tumor 
metabolic  demand,  hemoglobin  oxygen  saturation,  and  blood  flow 
would  provide  valuable  insight  into  the  mechanisms  and  distribution 
of  cycling  hypoxia  in  tumors.  Fluorescence  imaging  of  metabolic  de¬ 
mand  via  the  optical  redox  ratio  (fluorescence  intensity  of  FAD/ 
NADFH),  absorption  microscopy  of  hemoglobin  oxygen  saturation, 
and  Doppler  optical  coherence  tomography  of  vessel  morphology  and 
blood  flow  are  combined  to  noninvasively  monitor  changes  in  oxygen 
supply  and  demand  in  the  mouse  dorsal  skin  fold  window  chamber 
tumor  model  (human  squamous  cell  carcinoma)  every  6  h  for  36  h. 
Biomarkers  for  metabolic  demand,  blood  oxygenation,  and  blood 
flow  are  all  found  to  significantly  change  with  time  (p<0.05).  These 
variations  in  oxygen  supply  and  demand  are  superimposed  on  a  sig¬ 
nificant  (p<0.05)  decline  in  metabolic  demand  with  distance  from 
the  nearest  vessel  in  tumors  (this  gradient  was  not  observed  in  normal 
tissues).  Significant  (p<0.05),  but  weak  (r^0.5)  correlations  are 
found  between  the  hemoglobin  oxygen  saturation,  blood  flow,  and 
redox  ratio.  These  results  indicate  that  cycling  hypoxia  depends  on 
both  oxygen  supply  and  demand,  and  that  noninvasive  optical  imag¬ 
ing  could  be  a  valuable  tool  to  study  therapeutic  strategies  to  mitigate 
cycling  hypoxia,  thus  increasing  the  effectiveness  of  radiation  and 

chemotherapy.  ©  201 0  Society  of  Photo-Optical  Instrumentation  Engineers. 
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1  Introduction 

Hypoxic  tumors  are  more  aggressive  than  their  normoxic 
counterparts,  and  are  resistant  to  radiation  and 
chemotherapy,1"3  Thus,  a  major  goal  of  cancer  research  is  to 
gain  a  better  understanding  of  oxygen  supply  and  demand  in 
tumors.  Studies  of  oxygen  supply  and  demand  could  provide  a 
basis  for  improved  prognostic  and  treatment  approaches  for 
solid  tumors.  An  important  feature  of  tumor  hypoxia  is  its 
temporal  instability.  "Cycling  hypoxia”  has  been  described  as 
a  pattern  of  temporal  periodicity  between  hypoxic  and  reoxy¬ 
genated  states.  It  has  been  observed  in  many  tumor  types,4'6 
whereas  studies  measuring  p02  fluctuations  in  normal  tissue 
(muscle)  in  rats  and  mice  have  not  observed  significant 
fluctuations.7,8  One  primary  consequence  of  cycling  hypoxia 
is  upregulation  of  the  transcription  factor  hypoxia-inducible 
factor  1  (HI FI)  activity  to  a  level  that  supersedes  that  of 
chronic  hypoxia.  This  transcription  factor  controls  many  cel- 
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lular  functions  that  exacerbate  treatment  resistance  and  tumor 
aggressiveness.9 

Cycling  hypoxia  could  be  caused  by  changes  in  tumor  vas¬ 
cular  function  (oxygen  delivery)  and  metabolic  demand.  In 
tumors,  formation  of  new  vessels  and  remodeling  and  dropout 
of  existing  vessels  occurs  continuously,  resulting  in  temporal 
variations  in  blood  flow  and  blood  oxygenation.10  Oxygen 
consumption  rate  is  a  highly  dynamic  feature  of  oxygen  trans¬ 
port  in  tumors,  because  small  changes  in  demand  for  oxygen 
create  large  changes  in  the  extent  and  severity  of  hypoxia.11 
Longitudinal  imaging  of  tumor  metabolic  demand,  vessel 
morphology,  hemoglobin  oxygen  saturation,  and  blood  flow 
velocities  would  provide  valuable  insight  into  the  mechanisms 
and  distribution  of  cycling  hypoxia  in  tumors.  For  example, 
these  data  will  allow  for  the  study  of  tumor  response  to 
changing  hemodynamic  and  metabolic  conditions  induced  by 
experimental  and  traditional  cancer  therapies.  Such  informa¬ 
tion  can  also  be  used  to  test  and  refine  theoretical  models  for 
structural  adaptation  in  tumors,  which  have  successfully  pre¬ 
dicted  tumor  response  to  drug  and  radiation  therapy.12 
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Current  methods  for  measuring  oxygen  supply  on  the  mi¬ 
crovessel  level  include  microelectrode  measurements,7  and 
phosphorescence  lifetime  imaging  of  a  p02  calibrated  dye.13 
Microelectrodes  provide  a  sensitive  and  direct  measure  of 
p02,  but  are  inherently  invasive  and  limited  to  point  measure¬ 
ments.  Phosphorescence  lifetime  imaging  provides  informa¬ 
tion  on  the  2-D  distribution  of  oxygen  within  the  tumor,  but 
requires  the  injection  of  a  contrast  agent.  The  use  of  a  contrast 
agent  is  invasive  and  may  require  repeated  exposure  to  the 
dye  for  longitudinal  measurements,  which  further  complicates 
data  acquisition  and  interpretation.  Blood  flow  on  the  mi¬ 
crovessel  level  has  traditionally  been  measured  using  laser 
Doppler  flowmetry,7  the  dual-slit  method,13  or  video 
microscopy.14,15  Laser  Doppler  flowmetry  provides  a  relative 
(not  calibrated)  measurement  of  blood  flow  at  a  single  point, 
and  scanning  laser  Doppler  provides  an  average  flow  over  a 
tissue  area  without  accounting  for  contributions  from  indi¬ 
vidual  vessels.  The  dual-slit  technique  and  video  microscopy 
both  suffer  from  inaccurate  estimations  of  the  length  traveled 
by  red  blood  cells  because  the  length  of  a  3-D  vessel  is  pro¬ 
jected  onto  a  2-D  image,  resulting  in  imprecise  velocity  val¬ 
ues.  Video  microscopy  of  fluorescently  labeled  red  blood  cell 
flux  and  blood  flow  rate  can  be  done,  but  is  only  practical  for 
vessels  less  than  —30  ^m  in  diameter  because  individual  red 
blood  cells  are  difficult  to  visualize  in  larger  vessels  due  to 
absorption  of  fluorescent  light  by  hemoglobin. 

Previous  methods  for  measuring  oxygen  demand  in  tumors 
include  p02  microelectrode  measurements  coupled  with  the¬ 
oretical  simulations  of  oxygen  diffusion,16  an  isolated  tumor 
perfusion  system  combined  with  Fick's  principle 
calculations,1  and  cryospectrophotometric  microtechniques.18 
All  of  these  methods  are  invasive,  which  makes  longitudinal 
in  vivo  monitoring  difficult.  Positron  emission  tomography 
(PET)  of  fluoro-deoxyglucose  (FDG)  has  also  served  as  a 
valuable  marker  for  tumor  metabolism.19,20  Flowever,  PET  im¬ 
aging  is  low  resolution  and  expensive,  and  time  course  mea¬ 
surements  require  repeated  exposure  to  radioisotopes. 

Optical  imaging  of  tumor  oxygen  supply  and  demand  is 
attractive  because  these  methods  are  relatively  low  cost,  high 
resolution,  and  noninvasive,  if  performed  in  a  transparent 
window  chamber  preparation.  Optical  imaging  based  on  en¬ 
dogenous  tissue  contrast  does  not  require  exogenous  contrast 
agents,  and  is  thus  well  suited  for  long-term  monitoring  of 
cycling  hypoxia  in  tumors.  FI yperspectral  imaging  of  hemo¬ 
globin  oxygen  saturation21"23  coupled  with  3-D  Doppler  opti¬ 
cal  coherence  tomography  (OCT)  of  calibrated  microvessel 
blood  flow  rates24"27  allow  for  accurate,  dynamic,  and  nonin¬ 
vasive  imaging  of  oxygen  supply  in  tumors.28  Doppler  OCT 
also  allows  for  visualization  of  vessel  velocity  profiles  at  any 
vessel  cross  section  in  the  tumor,  thus  providing  information 
on  the  shear  rate  on  the  vessel  wall.29  Oxygen  demand  can  be 
optically  monitored  with  the  "redox  ratio,"  which  is  the  ratio 
of  the  inherent  fluorescence  intensity  of  the  metabolic  coen¬ 
zymes  FAD  and  N A D H  (the  primary  electron  acceptor  and 
donor,  respectively,  in  oxidative  phosphorylation.).30  This  op¬ 
tical  redox  ratio  provides  relative  changes  in  the  oxidation- 
reduction  state  in  the  cell.  The  redox  ratio  is  sensitive  to 
changes  in  the  cellular  metabolic  rate  and  vascular  oxygen 
supply.30-34 


The  goal  of  this  study  was  to  dynamically  image  changes 
in  oxygen  supply  and  demand  in  tumors  using  a  novel  com¬ 
bination  of  redox  imaging  (metabolic  demand)  and  hyper- 
spectral  imaging,  combined  with  Doppler  OCT  (oxygen  sup¬ 
ply).  The  results  reported  here  are  the  first  demonstration  of 
combined  noninvasive  imaging  of  oxygen  supply  and  meta¬ 
bolic  demand  in  tumors  in  vivo.  These  data  could  provide 
insight  into  the  sources  of  cycling  hypoxia  in  tumors,  and 
serve  as  the  basis  for  further  studies  that  will  incorporate 
therapeutic  strategies  to  mitigate  hypoxia  in  tumors  for  im¬ 
proved  treatment  response. 


2  Materials  and  Methods 

2.1  Multifunctional  Microscope 

Fluorescence,  hyperspectral,  and  OCT  imaging  were  con¬ 
ducted  on  the  same  microscope  to  noninvasively  monitor 
changes  in  metabolic  demand  and  hemodynamics  with  tumor 
growth.  All  images  were  collected  with  a  4x  objective  (NA 
=0.1).  FI  yperspectral  (hemoglobin  oxygen  saturation)  and 
SDOCT  (blood  flow  and  vascular  morphology)  arms  were 
connected  through  two  separate  base  ports  in  an  inverted  mi¬ 
croscope  (Carl  Zeiss  Axiovert  200).  2-D  hyperspectral  images 
were  collected  with  a  100-W  halogen  lamp  for  transillumina¬ 
tion,  and  detection  was  achieved  with  a  10-nm  bandwidth 
liquid  crystal  tunable  filter  (LCTF)  (CRI  Incorporated, 
Woburn,  Massachusetts)  placed  in  front  of  a  DVC  1412  CCD 
camera  (DVC  Company,  Austria,  Texas).23  Custom  software 
was  used  to  tune  the  filter  and  acquire  images  at  10-nm  in¬ 
crements  between  500  and  620  nm.  M  easurements  of  the  dark 
offset  and  transmission  through  a  neutral  density  filter  at  each 
wavelength  were  made  before  each  imaging  session.  Flemo- 
globin  saturation  images  were  calculated  by  applying  an  ex¬ 
tension  of  the  Beer-Lambert  law  to  the  wavelength-dependent 
attenuation  at  each  pixel  (assuming  oxygenated  and  deoxy- 
genated  hemoglobin  as  the  primary  absorbers,  along  with  a 
tissue  scattering  parameter),  and  then  solving  for  hemoglobin 
saturation  with  linear  least-squares  regression.23  The  system, 
software,  and  analysis  techniques  have  previously  been  vali¬ 
dated  on  liquid  phantoms  with  an  accuracy  of  approximately 
1%  and  in  vivoP 

2-D  fluorescence  imaging  was  conducted  in  the  epi- 
illumination  geometry  with  a  100-W  mercury  lamp.  N A D H 
images  were  collected  with  a  360-nm  bandpass  excitation  fil¬ 
ter  (40-nm  bandwidth)  and  a  390-nm  longpass  dichroic 
beamsplitter.  FAD  images  were  collected  with  a  470-nm 
bandpass  excitation  filter  (40-nm  bandwidth)  and  a  510-nm 
longpass  dichroic  beamsplitter.  Fluorescence  emission  was 
collected  with  the  same  LCTF  and  charge-coupled  device 
(CCD)  described  before  for  the  hyperspectral  arm.  The  LCTF 
was  set  to  480-nm  and  530-nm  emission  for  N  A  D  H  and  FAD 
imaging,  respectively.  Redox  images  were  calculated  by  di¬ 
viding  the  FAD  image  point  by  point  by  the  N A D H 
image.30,34  To  account  for  interanimal  variability  in  the  time- 
course  experiments,  the  redox  ratio  for  each  animal  was  nor¬ 
malized  the  first  time  point  of  the  first  region  of  interest 
within  each  animal.  Note  that  the  redox  image  represents  the 
redox  ratio  of  a  superficial  (—200  to  300  pm  thick)  tissue 
layer  probed  by  the  excitation-emission  wavelengths,35  and 
the  hemoglobin  saturation  images  show  vessels  from  a  simi- 
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larly  superficial  tissue  layer.22,23  OCT  volumes  indicate  that 
this  layer  contains  most  of  the  blood  vessels  in  the  field  of 
view.  The  spatial  resolution  of  hyperspectral  and  redox  imag¬ 
ing  is  defined  by  the  NA  of  the  objective  (5  [xm ). 

The  spectral  domain  OCT  arm  was  illuminated  by  a 
Ti:sapphire  laser  source  (Femtolasers,  Vienna,  Austria)  cen¬ 
tered  at  790  nm  with  a  90-nm  full  width  at  half-maximum. 
The  OCT  interferometer  was  common  path,  utilizing  the  top 
surface  of  the  window  chamber  cover  glass  as  the  reference 
reflection,  and  the  spectral  interferogram  was  detected  using  a 
custom-made  spectrometer  (spectral  resolution  and  coverage 
0.08  and  164  nm,  respectively)  with  a  2048  element  line  scan 
CCD  camera  (Atmel,  Aviiva,  San  Jose,  California).  The  OCT 
system  was  driven  by  software  that  controls  the  lateral  scan¬ 
ner  and  performs  data  acquisition  and  archiving  (Bioptigen, 
Incorporated,  Research  Triangle  park,  North  Carolina).  3-D 
blood  flow  images  were  collected  with  Doppler  OCT,  which 
measures  phase  changes  due  to  flowing  erythrocytes.24  Phase 
changes  were  calculated  from  multiple  A-scans  collected  at 
the  same  position  in  the  volume,  and  were  related  to  flow 
velocity  profiles  by  measuring  the  angle  of  incidence  in  the 
3-D  volume.  This  method  has  previously  been  validated  in 
flow  phantoms.36  Doppler  OCT  volumes  were  collected  over 
a  2x1.5  mm  area  with  250x  125  pixels  (with  a  resulting 
resolution  of  8  and  12  ^m  in  the  lateral  dimensions),  and 
from  2  mm  in  the  depth  dimension  with  1024  pixels,  with  a 
1-ms  integration  time  for  each  A-line.  Ten  repeated  A-scans  at 
each  location  were  collected  at  each  x-y  position  for  Doppler 
shift  estimation.  Cross  sectional  vessel  flow  profiles  were  fit 
to  a  second-order  polynomial  and  corrected  for  the  angle  of 
incidence  to  provide  velocity  in  millimeters  per  second.  Maxi¬ 
mum  blood  velocity  (Vmax)  was  determined  from  the  peak  of 
this  fit,  vessel  diameter  ( D )  was  determined  from  the  zero- 
crossings  of  the  fit,  and  the  shear  rate  on  the  vessel  wall  was 
determined  from  the  derivative  of  the  fit  (dvjdr)  evaluated  at 
the  vessel  wall,  assuming  a  Newtonian  fluid.29  Note  that  in 
vessels  with  diameters  30  to  60  p. m ,  the  velocity  profile  is 
expected  to  be  blunter  than  a  parabola  (this  reflects  the  non¬ 
continuum  behavior  of  blood  rather  than  non-Newtonian 
behavior).25,37-39  Thus,  the  estimated  wall  shear  rate  may  be 
an  underestimate  in  these  calculations.  However,  our  fit  re¬ 
sults  suggest  that  the  departure  from  a  parabola  is  small,  so 
the  second-order  polynomial  fit  is  a  reasonable  approxima¬ 
tion.  The  flow  in  the  vessel  is  the  average  velocity  in  the 
vessel  times  the  cross  sectional  area  of  the  vessel 
[Vmax/2*77*(£>/2)2]. 

2.2  In  Vivo  Imaging 

All  in  vivo  experiments  were  conducted  under  a  protocol  ap¬ 
proved  by  the  Duke  University  Institutional  Animal  Care  and 
Use  Committee.  A  titanium  window  chamber  was  surgically 
implanted  under  anesthesia  (ketamine  100  mg/kg  and  xyla- 
zine  10  mg/kg  IP)  on  the  back  of  athymic  nude  mice  (nu/nu, 
NCI,  Frederick,  Maryland).  A  window  chamber  tumor  was 
established  during  chamber  implantation  by  injecting  20  [xL 
of  a  single  cell  suspension  (8xl06  cells)  of  FaDu  (human 
squamous  cell  carcinoma)  cells  into  the  dorsal  skin  fold  prior 
to  placing  a  12-mm-diam  number  2  round  glass  coverslip 
(Erie  Scientific,  Portsmouth,  New  Hampshire)  over  the  ex¬ 
posed  tissue.  Animals  were  housed  in  an  on-site  housing  fa¬ 


cility  with  ad  libitum  access  to  food  and  water  and  standard 
12-h  light/dark  cycles. 

Redox  and  hemoglobin  saturation  imaging  were  validated 
on  mice  with  a  tumor-free  window  chamber.  Perturbation  ex¬ 
periments  were  performed  with  the  mice  under  ketamine/  xy- 
lazine  anesthesia  (maintained  at  body  temperature)  while 
breathing  air,  then  100%  oxygen  through  a  nose  cone,  and 
then  100%  nitrogen  through  a  nose  cone  at  1-L/min  flow 
rate. 

For  longitudinal  imaging  of  tumor  growth,  animals  were 
anesthetized  with  Isoflurane  (1.5%  with  oxygen)  and  main¬ 
tained  at  body  temperature.  Imaging  began  when  tumors 
reached  1  to  2  mm  in  diameter.  A  Doppler  OCT  volume 
(5-min  acquisition),  followed  by  hemoglobin  saturation  (30-s 
acquisition)  and  redox  (15-sec  acquisition)  imaging  were  col¬ 
lected  every  6  h  for  36  h.  Vessel  morphology  maps  were  used 
to  achieve  spatial  registration  between  each  imaging  session. 
Transition  between  microscope  modes  was  approximately 
1  min.  M  ice  were  allowed  to  recover  after  each  imaging  ses¬ 
sion. 

Note  that  no  corrections  were  made  to  the  data  to  account 
for  the  cardiac  cycle.  All  vessels  of  interest  were  microvessels 
(diameter  <70  /urn)  that  do  not  exhibit  large  changes  in  flow 
due  to  the  cardiac  cycle.  Previous  studies  have  used  Fourier 
analysis  of  blood  flow  and  p02  in  rats  to  show  low-frequency 
fluctuations  (<2  cycles/min)  in  both  tumor  and  muscle  tis¬ 
sues,  which  is  much  lower  than  the  measured  respiration  rate 
and  heart  rate  (~50  to  60  and  300  to  400  cycles/min, 
respectively.7 

2.3  Statistical  Analysis 

Doppler  OCT  and  hyperspectral  data  were  collected  over  the 
36-h  time  course  from  a  total  of  five  mice,  and  redox  images 
were  also  collected  over  the  time  course  from  three  of  these 
mice.  Three  to  six  regions  of  interest  were  chosen  within  each 
tumor  based  on  visual  inspection  of  the  hemoglobin  saturation 
time  course  images.  Vessel  regions  that  demonstrated  fluctua¬ 
tions  in  hemoglobin  saturation  over  the  time  course,  and  that 
were  spatially  distributed  across  the  surface  of  the  tumor, 
were  chosen  as  regions  of  interest.  The  mean  and  standard 
deviation  of  the  hemoglobin  saturation  were  measured  at 
these  vessel  regions  of  interest,  over  an  area  within  the  inner 
diameter  of  the  vessel  and  a  slab  thickness  of  no  more  than 
20-  ytim.  The  Doppler  velocity  profile  was  measured  at  the 
same  vessel  cross  section  of  interest.  The  mean  and  standard 
deviation  of  the  redox  ratio  was  measured  from  a 
25  to  60-^t.m-diam  circular  area  adjacent  to  the  vessel.  The 
redox  area  of  interest  was  chosen  to  exclude  blood  vessels  and 
to  remain  less  than  100  /x m  in  diameter  (the  length  of  oxygen 
diffusion  in  tissue.40  A  total  of  21  regions  of  interest  were 
analyzed  over  the  time  course,  with  14  regions  of  interest  that 
include  redox  ratio  measurements.  One  of  the  animals  for 
which  redox  images  were  collected  included  only  six  time 
points  for  the  four  regions  of  interest  within  that  animal,  due 
to  tumor  hemorrhage  at  the  last  time  point. 

Pearson's  correlation  coefficients  were  used  to  determine 
linear  correlations  between  blood  oxygenation  (hemoglobin 
saturation),  blood  flow  (maximum  velocity,  flow,  vessel  inner 
diameter,  and  shear  rate),  and  metabolic  demand  (redox  ratio). 
Significant  correlations  (/? < 0.05)  are  reported  for  the  pooled 
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data  (all  regions  of  interest  and  time  points  grouped  together), 
for  which  there  were  n  = 143  samples  for  the  blood  oxygen¬ 
ation  and  blood  flow  comparisons,  and  n= 94  samples  for  the 
metabolic  demand  comparisons.  As  described  before,  the  re¬ 
dox  ratio  was  normalized  to  the  first  time  point  of  the  first 
region  of  interest  within  each  animal.  Thus,  blood  oxygen¬ 
ation  and  blood  flow  biomarkers  were  also  normalized  to  the 
first  time  point  of  the  first  region  of  interest  within  each  ani¬ 
mal  for  comparisons  with  the  redox  ratio  only. 

Next,  a  linear  mixed  model  was  used  to  analyze  repeated 
measures.  This  approach  separates  the  interanimal  variation 
from  the  total  variation  (thus  increasing  statistical  power)  to 
determine  whether  the  biomarkers  (i.e„  hemoglobin  satura¬ 
tion,  redox  ratio,  maximum  velocity,  vessel  inner  diameter, 
flow,  and  shear  rate)  change  over  time.  Let  Yijk  denote  the 
value  of  a  biomarker  from  the  ;'th  animal  in  the / th  region  of 
interest  (ROI)  at  the  &’th  time  point.  Assume  Yijk  follows  the 
linear  mixed  model  Yijk=/j.+tk+bij+eijk  (model  1);  where 
fi+tk  is  the  time  effect,  is  the  animal  specific  random 
effect,  and  sijk  are  intra-animal  errors.  The  time  effect  was 
tested  by  the  F-statistic.  A  second  mixed  model  was  used  to 
address  the  association  of  two  variables  adjusting  for  the  co¬ 
variate  time.  Model  2  ( Yljk = +Xljk+tk Ybij Jrsl]k)  includes 
the  covariate  X,  and  the  linear  association  of  biomarker  X  on 
the  behavior  of  another  biomarker  Y  was  tested  with  the 
F-statistic.  Associations  were  tested  between  blood  oxygen¬ 
ation  (hemoglobin  saturation),  blood  flow  (maximum  velocity, 
flow,  vessel  inner  diameter,  and  shear  rate),  and  metabolic 
demand  (redox  ratio)  biomarkers. 

3  Results 

The  results  of  perturbation  experiments  in  a  normal  mouse 
window  chamber  (no  tumor  cells  were  injected)  are  shown  in 
Fig.  1.  To  segment  the  vessels  from  background,  the  hemo¬ 
globin  saturation  images  were  thresholded  by  the  total  absorp¬ 
tion  image.  The  same  mask  was  negatively  applied  to  the 
redox  images.  Note  that  there  are  two  distinct  peaks  in  the 
hemoglobin  saturation  histogram  for  air  and  oxygen  breathing 
[Fig.  1(b)],  due  to  arteriolar  and  venous  blood  oxygenation. 
Nitrogen  and  oxygen  breathing  resulted  in  a  decrease  and  in¬ 
crease,  respectively,  in  both  the  redox  ratio  and  the  hemoglo¬ 
bin  saturation  compared  to  air  breathing,  as  expected.  The 
results  in  Fig.  1  were  confirmed  in  a  second  mouse,  and  are 
consistent  with  previous  reports.23,33  These  perturbation  ex¬ 
periments  indicate  that  NADH  and  FAD  signals  have  been 
properly  isolated  for  redox  imaging. 

The  hemoglobin  saturation  [Fig.  1(a)]  and  redox  [Fig. 
1(b)]  images  collected  from  one  mouse  over  the  36-h  time 
course  are  shown  in  Fig.  2.  Longitudinal  oxygen  gradients  are 
apparent  in  the  hemoglobin  saturation  images,  with  decreased 
hemoglobin  saturation  levels  within  the  tumor  [Fig.  2(a),  at 
the  1]  compared  with  the  tumor  periphery  [Fig.  2(a),  at  the  2], 
These  longitudinal  gradients  have  been  previously  observed 
in  studies  employing  phosphorescence  lifetime  imaging  of 
/?02.41,42  The  redox  images  indicate  that  the  redox  ratio  within 
the  tumor  is  higher  in  areas  with  large-diameter  oxygenated 
vessels  [Fig.  2(b),  at  the  3],  Both  the  hemoglobin  saturation 
and  the  redox  ratio  appear  to  vary  with  time. 

Figure  3  shows  an  en-face  view  of  the  3-D  Doppler  OCT 
volume  (manually  segmented  in  Amira  software,  Mercury 


(a) 


Fig.  1  Hemoglobin  saturation  (HbSat)  and  redox  images  taken  from  a 
nontumor-bearing  window  chamber  of  a  mouse  breathing  oxygen,  air, 
and  nitrogen  (a)  at  1  L/min  through  a  nose  cone.  To  segment  the 
vessels  from  background,  the  hemoglobin  saturation  images  were 
thresholded  by  the  total  absorption  image.  The  same  mask  was  nega¬ 
tively  applied  to  the  redox  images.  Histograms  of  (b)  the  hemoglobin 
saturation  and  (c)  redox  images  show  that  the  hemoglobin  saturation 
and  redox  ratio  increase  with  oxygen  breathing  and  decrease  with 
nitrogen  breathing,  as  expected.  Image  sizes  are  2X  1.5  mm. 


Systems,  Chelmsford,  M  assachusetts)  from  the  zero-hour  time 
point  [Fig.  3(a)]  of  the  same  animal  shown  in  Fig.  2.  Plots  of 
the  cross  sectional  vessel  velocity  profile  from  one  vessel  [ar¬ 
row,  Fig.  3(a)]  are  shown  along  with  the  second-order  poly¬ 
nomial  fit  and  the  R-squared  value  of  the  fit  for  each  time 
point  [Figs.  3(b)- 3(h)].  The  shape  and  peak  of  the  velocity 
profile  changes  with  time,  indicating  changes  in  vessel  diam¬ 
eter,  shear  rate,  and  maximum  vessel  velocity  over  the  36-h 
time  course. 

Figure  4  shows  quantitative  measures  of  the  redox  ratio 
[Fig.  4(b)],  hemoglobin  saturation  [Fig.  4(c)],  vessel  maxi¬ 
mum  velocity  [Fig.  4(d)],  vessel  inner  diameter  [Fig.  4(e)], 
flow  [Fig.  4(f)],  and  shear  rate  [Fig.  4(g)]  over  the  36-h  time 
course  in  six  regions  of  interest  [Fig.  4(a)]  from  the  same 
animal  shown  in  Figs.  2  and  3.  These  plots  indicate  that  all 
variables  change  with  time  within  the  same  region  of  interest 
within  the  same  animal.  All  regions  of  interest  appear  to  vary 
together  with  time  for  the  redox  ratio  [Fig.  4(b)]  and  for  the 
hemoglobin  saturation  [Fig.  4(c)].  However,  each  region  of 
interest  appears  to  vary  differently  with  time  compared  with 
the  other  five  regions  of  interest  for  the  blood  flow  parameters 
[Figs.  4(d)- 4(g)], 

The  spatial  distribution  of  the  redox  ratio  with  respect  to 
the  nearest  vessel  is  shown  in  Fig.  5  for  tumors  and  normal 
tissues.  There  is  a  significant  (/? < 0.05)  decrease  in  the  redox 
ratio  for  points  greater  than  30  /xm  from  the  nearest  vessel, 
compared  to  the  point  closest  to  the  nearest  vessel  (point 
zero).  This  is  in  agreement  with  previous  reports  that  mea¬ 
sured  p02  as  a  function  of  distance  from  the  nearest 
vessel.15,16  However,  the  normal  tissues  showed  no  significant 
increase  or  decrease  in  the  redox  ratio  as  a  function  of  dis¬ 
tance  from  the  nearest  vessel. 
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Fig.  2  (a)  Hemoglobin  saturation  and  (b)  redox  images  taken  from 
one  mouse  over  the  36-h  time  course.  The  time  of  imaging  in  hours  is 
indicated  in  the  lower  left  corner  of  each  hemoglobin  saturation  im¬ 
age.  For  the  purposes  of  this  figure  only,  the  hemoglobin  saturation 
images  were  thresholded  by  the  total  absorption  image,  and  then  by  a 
1 5-/xm-diam  circular  kernel.  The  kernel  excluded  pixels  with  an  ab¬ 
sorption  value  smaller  than  the  average  of  the  15-/xm-diam  kernel. 
The  mask  was  applied  to  hemoglobin  saturation  images  to  isolate 
vessels  from  the  nonvascular  background,  and  the  same  vessel  mask 
was  negatively  applied  to  the  redox  images.  The  tumor  boundary  is 
indicated  by  a  white  dashed  line  in  (b).  Image  sizes  are  2  X  1 .5  mm. 


The  relationship  between  the  biomarkers  for  blood  oxy¬ 
genation  (hemoglobin  saturation),  blood  flow  (maximum  ve¬ 
locity,  flow,  vessel  inner  diameter,  and  shear  rate),  and  meta¬ 
bolic  demand  (redox  ratio)  were  tested  with  the  Pearson's 
correlation  coefficient  (Table  1)  for  the  pooled  data  (all  re¬ 
gions  of  interest  and  time  points  from  all  animals  grouped 
together).  The  correlation  coefficients  (r)  in  Table  1  are 
shown  for  biomarker  relationships  for  which  a  significant  (p 
<0.05)  linear  correlation  was  found.  The  hemoglobin  satura¬ 
tion  is  positively  correlated  with  the  maximum  velocity,  flow, 
and  shear  rate  (HbSat  line  of  Table  1),  and  the  redox  ratio  is 
positively  correlated  with  the  hemoglobin  saturation,  maxi¬ 
mum  velocity,  flow,  and  diameter.  All  statistically  significant 
correlations  are  weak  (r«0.5),  which  indicates  that  many 
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Fig.  3  An  en-face  view  of  the  3-D  Doppler  OCT  volume  (manually 
segmented  in  Amira  software,  Mercury  Systems)  from  (a)  the  zero- 
hour  time  point  of  the  same  animal  shown  in  Fig.  2  (2Xl.5-mm 
imaged  area,  0.6  mm  depth).  Plots  of  the  cross  sectional  vessel  veloc¬ 
ity  profile  from  one  vessel  [arrow  in  (a)]  are  shown  along  with  the 
second-order  polynomial  fit  and  the  R-squared  value  of  the  fit  for  each 
time  point  (b)  through  (h).  The  time  of  imaging  in  hours  is  indicated  in 
the  upper  right  corner  of  each  plot. 


factors  contribute  to  the  variability  of  each  biomarker,  as 
expected.11  The  linear  mixed  model  determined  that  all  biom¬ 
arkers  significantly  change  (p<0.05)  with  time  (Fig.  6),  and 
that  the  vessel  diameter  is  associated  with  the  redox  ratio  of 
the  adjacent  tissue  (p < 0.05)  after  adjusting  for  the  effect  of 
time.  No  other  biomarker  pairs  were  found  to  have  significant 
associations  (p >  0.05)  when  adjusting  for  the  effect  of  time. 

4  Discussion 

Multifunctional  optical  imaging  of  oxygen  supply  and  de¬ 
mand  has  several  advantages  over  traditional  techniques.  The 
combination  of  OCT,  absorption,  and  fluorescence  imaging 
provides  quantitative  information  on  blood  flow,  blood  oxy¬ 
genation,  and  metabolic  demand,  respectively,  in  living  tis¬ 
sues.  This  combination  of  optical  imaging  modalities  is  non- 
invasive,  relatively  inexpensive,  and  fast.  These  attributes  are 
well  suited  for  long-term  repeated  measures  in  vivo.  We  have 
demonstrated  the  attributes  of  this  imaging  platform  for  the 
study  of  cycling  hypoxia  in  tumors  in  the  dorsal  skin  fold 
mouse  window  chamber.  This  is  the  first  demonstration  of 
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Fig.  4  Quantitative  measures  of  (b)  the  redox  ratio,  (c)  hemoglobin 
saturation,  (d)  vessel  maximum  velocity,  (e)  vessel  inner  diameter,  (f) 
flow,  and  (g)  shear  rate  over  the  36-h  time  course  in  (a)  six  regions  of 
interest  from  the  same  animal  shown  in  Figs.  2  and  3.  The  mean  and 
standard  error  are  plotted  for  (b)  the  redox  ratio  and  (c)  hemoglobin 
saturation,  and  the  mean  (d)  through  (g)  is  plotted  for  the  remaining 
variables. 


combined  noninvasive  monitoring  of  oxygen  supply  and 
metabolic  demand  in  tumors  in  vivo. 

All  biomarkers  (hemoglobin  saturation,  redox  ratio,  maxi¬ 
mum  velocity,  vessel  inner  diameter,  flow,  and  shear  rate) 
significantly  changed  with  time  (/? < 0.05)  over  the  36-h  im¬ 
aging  period,  which  supports  previous  reports  that  blood  oxy¬ 
genation,  blood  flow,  and  metabolic  demand  all  contribute  to 
cycling  hypoxia  in  tumors.3  These  variations  in  oxygen  sup¬ 
ply  and  demand  are  superimposed  on  a  decline  in  metabolic 
demand  with  distance  from  the  nearest  vessel  in  tumors  (Fig. 
5),  in  agreement  with  previous  interstitial  p02  microelectrode 
measurements.15,16  Our  results  indicate  a  positive  correlation 
between  the  blood  flow,  blood  oxygenation,  and  metabolic 
demand  in  tumors  (Table  1),  but  the  correlation  between  these 
biomarkers  is  weak  and  suggests  that  no  single  biomarker  can 
fully  explain  the  behavior  of  another.11  W  hen  adjusting  for  the 
effect  of  time,  the  vessel  diameter  was  found  to  be  associated 
with  the  redox  ratio  of  the  adjacent  tissue,  which  indicates 
that  vessel  diameter  plays  a  key  role  in  tumor  metabolic  ac¬ 
tivity  and  cycling  hypoxia  in  tumors.43  Future  studies  will 
exploit  these  biomarkers  to  investigate  functional  changes  in 


Fig.  5  Normalized  redox  ratio  as  a  function  of  distance  from  the  near¬ 
est  vessel  for  tumors  (compiled  from  n  =  94  line  profiles  corresponding 
to  the  94  regions  of  interest  from  the  three  redox  animals  across  the 
entire  tumor  time  course),  and  for  normal  tissues  (compiled  from  n 
=  36  line  profiles  from  six  regions  of  interest  from  each  of  six  normal 
animals).  All  profiles  were  normalized  to  the  point  closest  to  the  ves¬ 
sel  before  averaging  across  all  line  profiles.  For  the  tumors,  there  is  a 
significant  (p<0.05)  decrease  in  the  redox  ratio  between  the  point 
closest  to  the  vessel  (point  zero)  and  all  points  greater  than  30  pm 
from  point  zero.  Normal  tissues  did  not  show  a  significant  increase  or 
decrease  in  the  redox  ratio  as  a  function  of  the  distance  from  the 
nearest  vessel.  Error  bars  are  standard  error. 


oxygen  supply  and  demand  in  response  to  therapeutic  strate¬ 
gies.  For  example,  it  is  possible  that  antiangiogenic  therapies 
may  stabilize  cycling  hypoxia  by  reducing  the  magnitude  of 
changes  in  vascular  diameter  and  blood  flow  rate  over  time.3 
It  is  anticipated,  however,  that  antiangiogenic  therapies  would 
have  less  influence  on  redox  ratio.  The  combination  of  tech¬ 
nologies  presented  here  would  permit  a  direct  test  of  this  hy¬ 
pothesis. 

Extracted  biomarkers  including  the  hemoglobin 
saturation,44,45  maximum  velocity  and  vessel  diameter,14  and 
shear  rate,46  are  all  in  qualitative  agreement  with  previous 
studies  that  employed  different  techniques  in  window  cham¬ 
ber  tumors,  as  well  as  a  previous  study  that  employed  the 
same  techniques  in  a  different  window  chamber  tumor 
model.28  Variations  in  red  blood  cell  flux  and  p02,15  hemo¬ 
globin  oxygen  saturation,23  and  diameter44  have  previously 


Table  1  Pearson  correlation  coefficients  for  hemoglobin  saturation 
(FHbSat)  versus  maximum  velocity  (V^),  flow,  diameter,  and  shear 
rate  (top  line)  and  for  the  normalized  redox  ratio,  normalized  hemo¬ 
globin  saturation,  normalized  maximum  velocity,  normalized  flow, 
normalized  diameter  and  normalized  shear  rate  (bottom  line).  The 
data  are  compiled  from  all  regions  of  interest  and  time  points  grouped 
together,  with  n=143  observations  for  the  top  line  and  n  =  94  obser¬ 
vations  for  the  bottom  line.  The  correlation  coefficient  is  listed  for 
comparisons  with  p<0.05. 


v 

r  max 

(mm/s) 

Flow 
(mm3/ s) 

Diameter 

(pm) 

Shear  rate 

(1/s) 

HbSat  (%) 

0.49 

0.37 

p>0.05 

0.49 

HbSat 

y  max 

Flow 

Diameter 

Shear  rate 

Redox 

0.50 

0.28 

0.32 

0.43 

p>0.05 
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Fig.  6  Longitudinal  variations  in  all  biomarkers  as  a  function  of  time. 
Values  at  each  time  point  were  averaged  across  all  regions  of  interest 
and  animals  for  the  hemoglobin  saturation  (HbSat),  flow,  diameter, 
maximum  velocity  (Vmax),  and  shear  rate  (n= 21),  as  well  as  for  the 
redox  ratio  (n=  14).  Note  that  the  36-h  time  point  has  four  fewer 
samples  (see  Sec.  2).  The  linear  mixed  model  determined  that  all 
biomarkers  significantly  change  (p<  0.05)  with  time,  and  that  the  ves¬ 
sel  diameter  is  associated  with  the  redox  ratio  of  the  adjacent  tissue 
(p<0.05)  after  adjusting  for  the  effect  of  time.  Error  bars  are  standard 
error. 


been  observed  in  window  chamber  tumors  over  time  scales 
ranging  from  one  hour  to  several  days. 

Previous  studies  have  employed  different  methods  to  in¬ 
vestigate  the  relationship  between  some  of  the  biomarkers 
presented  in  this  work.  Intravascular p02  and  blood  flow  rate, 
as  well  as  interstitial  p02  and  red  blood  cell  velocity  or  vessel 
diameter,  were  found  to  have  no  correlation  in  a  previous 
study.13  Helmlinger  et  al.  used  phosphorescence  lifetime  im¬ 
aging  to  quantify  p02  and  video  microscopy  to  quantify  red 
blood  cell  velocity  and  vessel  diameter.  Although  hemoglobin 
oxygen  saturation  and  p02  are  related,  they  are  not  the  same 
endpoint,  and  video  microscopy  provides  a  2-D  projection 
rather  than  an  absolute  measure  of  blood  velocity.  The  study 
design  and  tumor  models  employed  in  the  previous  and  cur¬ 
rent  studies  also  differ.  Previous  studies  employing  the  same 
techniques  as  the  current  study  have  found  a  significant  cor¬ 
relation  between  the  hemoglobin  saturation  and  blood  veloc¬ 
ity,  blood  flow,  and  shear  rate  in  a  different  tumor  model.28 
Additional  studies15  indicate  that  vessel  remodeling  with  tu¬ 
mor  growth  alters  flow  distributions,  and  the  relationship  be¬ 
tween  blood  flow  and  tissue  oxygenation  in  tumors  is  hetero¬ 
geneous. 

Demand  for  oxygen  in  tumor  cells  is  created  by  the  normal 
respirative  process  of  the  cells.32  The  redox  ratio  is  a  measure 
of  cell  respiration  (oxygen  demand),30  while  the  hemoglobin 
oxygen  saturation  is  one  component  of  oxygen  supply.  There¬ 
fore,  it  is  intuitive  that  the  hemoglobin  saturation  of  a  vessel  is 
somewhat  related  to  the  redox  ratio  of  the  surrounding  tissue. 
However,  we  do  not  expect  a  one-to-one  correlation  between 
the  redox  ratio  and  hemoglobin  saturation.  Oxygen  delivery  is 
affected  by  the  hemoglobin  saturation  as  well  as  the  vascular 
architecture,  hematocrit,  and  blood  flow.32  The  redox  ratio  is 
indicative  of  changes  in  oxygen  supply  as  well  as  metabolite 
availability  and  mitochondrial  membrane  potential.47  Thus, 
the  redox  ratio,  blood  flow  and  hemoglobin  saturation  provide 
complementary  information. 

Figure  5  shows  a  decrease  in  the  redox  ratio  with  increas¬ 
ing  distance  from  the  nearest  vessel  in  tumors,  and  no  such 
decrease  in  the  redox  ratio  with  distance  from  the  nearest 
vessel  in  normal  tissues.  We  hypothesize  that  this  tumor  gra¬ 
dient  is  due  to  increased  metabolic  activity  coupled  with  com¬ 


promised  microvascular  oxygen  transport  in  tumors  compared 
to  normal  tissues.  Future  studies  will  investigate  this  hypoth¬ 
esis. 

In  this  study,  absorption  and  fluorescence  microscopy  were 
combined  with  OCT  to  unravel  the  contributions  of  blood 
flow,  blood  oxygenation,  and  metabolic  demand  on  cycling 
hypoxia.  These  findings  could  guide  the  design  and  develop¬ 
ment  of  therapeutic  strategies  to  mitigate  cycling  hypoxia,  and 
thus  reduce  tumor  resistance  to  radiation  and  chemotherapies. 
Second,  this  study  experimentally  demonstrates  the  attributes 
of  our  multifunctional  optical  imaging  approach  for  mapping 
oxygen  supply  and  demand  in  vivo.  This  approach  has  broad 
applicability  in  the  study  of  other  diseases  that  are  affected  by 
oxygen  supply  and  demand,  such  as  cerebral  hypoxia  (stroke), 
cardiac  ischemia,  and  Alzheimer's  disease.  Moreover,  these 
processes  can  be  studied  with  multifunctional  optical  imaging 
in  vivo  and  do  not  require  contrast-enhancing  dyes. 
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Abstract.  The  influence  of  the  tumor  microenvironment  and  hypoxia 
plays  a  significant  role  in  determining  cancer  progression,  treatment 
response,  and  treatment  resistance.  That  the  tumor  microenvironment 
is  highly  heterogeneous  with  significant  intratumor  and  intertumor  vari¬ 
ability  presents  a  significant  challenge  in  developing  effective  cancer 
therapies.  Critical  to  understanding  the  role  of  the  tumor  microenvi¬ 
ronment  is  the  ability  to  dynamically  quantify  oxygen  levels  in  the 
vasculature  and  tissue  in  order  to  elucidate  the  roles  of  oxygen  supply 
and  consumption,  spatially  and  temporally.  To  this  end,  we  describe 
the  use  of  hyperspectral  imaging  to  characterize  hemoglobin  absorp¬ 
tion  to  quantify  hemoglobin  content  and  oxygen  saturation,  as  well  as 
dual  emissive  fluorescent/phosphorescent  boron  nanoparticles,  which 
serve  as  ratiometric  indicators  of  tissue  oxygen  tension.  Applying  these 
techniques  to  a  window-chamber  tumor  model  illustrates  the  role  of 
fluctuations  in  hemoglobin  saturation  in  driving  changes  in  tissue  oxy¬ 
genation,  the  two  being  significantly  correlated  (r  =  0.77).  Finally,  a 
green-fluorescence-protein  reporter  for  hypoxia  inducible  factor-1  (HI F- 
1 )  provides  an  endpoint  for  hypoxic  stress  in  the  tumor,  which  is  used  to 
demonstrate  a  significant  association  between  tumor  hypoxia  dynamics 
and  HIF-1  activity  in  an  in  vivo  demonstration  of  the  technique.  ©2010 
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Tumor  hypoxia  plays  an  important  role  in  cancer  progression, 
metastasis,  and  the  effectiveness  of  therapies.1-3  It  is  a  highly 
dynamic  process  and  the  temporal  dynamics  can  play  a  signifi¬ 
cant  role  in  treatment  resistance  and  tumor  aggressiveness.2,4-10 
A  recent  trend  in  therapeutic  strategies  is  to  target  tumor  an¬ 
giogenesis  in  order  to  modify  the  microenvironment  and  make 
tumors  more  sensitive  to  conventional  therapies.11  The  ability 
to  measure  pC>2  and  image  tumor  hypoxia  on  the  microscopic 
level  is  essential  to  further  understanding  microenvironmen¬ 
tal  factors  and  modifying  therapies,  depending  on  treatment 
response.  Characterizing  the  local  microenvironment,  partic¬ 
ularly  hypoxia,  longitudinally,  on  a  microscopic  scale  would 
be  ideal.  This  is  a  challenging  problem,  however,  and  one  for 
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which  no  established  imaging  technique  exists.  Many  methods 
have  been  developed  that  range  from  histological  staining  with 
reducible  nitroimidazole  reagents,  protein  markers,  and  green 
fluorescence  protein  (GFP)  or  lucif erase  reporter  genes,  to  elec¬ 
trochemical  probe  methods,  and  optical,  nuclear  medicine,  and 
magnetic  resonance  techniques12  although  some  methods  are 
only  suitable  for  the  laboratory,  others  have  been  translated  into 
the  clinic,  and  all  have  their  advantages  and  drawbacks.13,14 
The  key  point  is  that  there  is  no  easily  implemented  technique 
for  imaging  vascular  oxygenation  (hemoglobin  oxygenation)  in 
combination  with  tissue  oxygen  tension  on  the  microcirculatory 
level.  We  present  an  approach  based  on  hyperspectral  imaging  of 
hemoglobin  absorption  and  dual  emissive  boron  nanoparticles 
(BNPs).  The  combination  of  these  techniques  enables  quantifi¬ 
cation  of  both  oxygen  supply  (hemoglobin  oxygen  saturation) 
and  tissue  oxygen  tension  (BNPs  ratiometric  sensing). 
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Recently,  we  reported  a  new  kind  of  light-emitting  biomate¬ 
rial  that  offers  many  advantages  for  optical  hypoxia  imaging15 
Iodide- substituted  difluoroboron  dibenzoylmethane-polylactide 
(BF2dbm(I)PLA)  is  multiemissive,  exhibiting  both  short-lived 
fluorescence  (F)  and  long-lived  phosphorescence  (P)  after  1- 
or  2-photon  excitation.  Especially  unique  is  the  presence  of 
P  at  room  temperature  (and  37°C)  in  these  PL  A  biomateri¬ 
als,  even  as  BNPs  in  aqueous  environments,16  which  makes 
them  useful  for  O2  sensing  via  a  dynamic  quenching  mecha¬ 
nism.  Many  pC>2  sensors  are  based  on  lifetime  (r),  because  r 
is  a  concentration-independent  variable,  but  its  measurement 
requires  specialized  instrumentation  not  typically  present  in 
routine  fluorescence  detection  systems.  More  straightforward, 
however,  is  ratiometric  sensing.  A  simple  ratio  of  F/P  inten¬ 
sities  at  the  respective  emission  maxima  provides  information 
about  relative  O2  levels,  whereas  calibration  allows  for  abso¬ 
lute  pC>2  measurement.  Ratiometric  methods  also  eliminate  the 
effects  of  concentration  and  fluctuations  in  light  intensity  or 
detector  sensitivity.  Previously,  ratiometry  had  been  achieved 
with  ternary  mixtures  of  a  cyanine  dye  standard  and  an  O2- 
sensitive  Pt  porphyrin  phosphor  in  a  sol  gel  matrix.17  But  single¬ 
component,  dual-emissive  BF2dbm(I)PLA  offers  the  additional 
advantages  of  greater  sample  homogeneity,  minimal  dye  leach¬ 
ing,  biocompatibility,  and  simplified  processing;  because  the 
nanoparticle  serves  as  both  the  standard  (F)  and  the  oxygen 
sensor  (P)  at  once,  should  a  sensor  molecule  degrade,  intrinsi¬ 
cally  coupled  F  and  P  are  equally  affected.  Furthermore,  long 
P  lifetimes  (r  =  ms;  compare  to  F:  r  =0.5  ns)  corre¬ 
late  with  high  oxygen  sensitivity,  making  these  materials  es¬ 
pecially  well  suited  for  hypoxia  imaging.  The  second  parame¬ 
ter  of  interest  is  vascular  oxygenation,  which  can  be  obtained 
using  hyperspectral  imaging  of  transmitted  light  through  the 
window  chamber  to  characterize  hemoglobin  absorption  as  a 
function  of  wavelength.  From  these  data,  it  is  possible  to  ex¬ 
tract  the  hemoglobin  oxygen  saturation  within  the  microvas¬ 
culature,  as  well  as  the  product  of  hemoglobin  concentration 
and  path  length,  which  is  proportional  to  the  total  attenuation 
through  the  tissue.  Finally,  a  HIF-1  reporter  construct  is  ex¬ 
pressed  in  the  4T1  cell  line  used  in  this  study,  which  provides  a 
direct  measure  of  HIF- 1  activity  through  induced  expression  of 
GFP.18  This  enables  quantification  of  the  biologic  consequences 
of  hypoxia  on  an  important  hypoxia-associated  pathway  that 
has  been  shown  to  influence  treatment  resistance  and  tumor 
aggressiveness.10 


1  Methods  and  Materials 

1 .1  BNP  Preparation 

The  dye-polymer  conjugate15  BF2dbm(I)PLA  and  boron 
nanoparticle  aqueous  suspensions16  were  fabricated  accordingly 
to  previously  reported  methods  at  a  concentration  of  1  mg/ml. 
Specifically,  the  polymer  was  generated  by  solvent-free  lactide 
ring  opening  polymerization  using  the  primary  alcohol  func¬ 
tionalized  BF2dbm(I)OH  initiator  and  tin  octoate  as  the  cat¬ 
alyst.  Nanoprecipitation  of  N,N-dimethylformamide  solutions 
of  BF2dbm(I)PLA  into  water,  followed  by  dialysis  and  pas¬ 
sage  through  filter  paper,  yields  nanoparticles  ~ 80- 100  nm 
diam. 


1 .2  Microscope 

A  Zeiss  MPS  intravital  microscope  was  used  for  ah  imaging, 
using  a  2.5X  objective.  A  DAPI  excitation  filter  cube  was  used 
for  excitation,  and  a  liquid-crystal  tunable  filter  was  used  to 
isolate  the  emission  wavelengths  (Varispec,  Cri  Inc.,  Woburn, 
MA).  This  filter  has  a  10-nm  bandwidth  with  an  electronically 
controllable  center  wavelength.  A  CCD  camera  (AM  14 12,  Dig¬ 
ital  Video  Camera  Company,  Austin,  TX)  was  then  used  to 
record  images.  For  acquisition  of  BNP  data,  the  DAPI  excita¬ 
tion  filter  and  mercury  lamp  source  (HBOIOO,  Carl  Zeiss  AG, 
Oberkochen,  Germany)  were  used  to  excite  the  sample,  with 
the  emission  recorded  from  425  to  575  nm  in  5-nm  increments 
using  1-s  exposures.  This  resulted  in  an  31-s  acquisition  time 
for  oxygen  tension  maps.  These  were  acquired  in  increments 
of  7.5  min,  with  brightfield  hemoglobin  saturation  acquisitions 
acquired  continuously  in-between  these  (GFC  Mass  Flow  Con¬ 
troller,  Aalborg  Inc.,  Orangeburg,  NY).  This  was  done  to  mini¬ 
mize  exposure  of  tissue  to  UV  light  during  the  hour  acquisition 
time.  For  brightfield  imaging  used  for  calculation  of  hemoglobin 
saturation,  an  unfiltered  halogen  white-light  source  was  used  in 
transmission  mode  through  the  window,  with  emission  recorded 
from  520  to  620  nm  in  10-nm  increments. 

1 .3  BNP  calibration 

The  BNP  solution  was  imaged  in  a  1-cm  quartz  cuvette  with 
continuous  bubbling  of  gas.  The  gas  was  a  mixture  of  nitrogen 
and  oxygen,  with  flow  rates  controlled  by  mass  flow  controllers 
(Aalborg,  Inc.).  The  solution  was  initially  purged  of  oxygen  by 
bubbling  with  100  N2  for  30  min  and  then  imaged  using  the 
setting  described  above.  The  oxygen  flow  was  then  incremen¬ 
tally  increased  to  provide  a  range  of  oxygen  concentration  from 
0  to  20%,  with  images  taken  at  each  increment  after  a  5 -min 
stabilization  period. 

Calibration  of  BNP  image  data  to  absolute  oxygen  tension 
was  performed  as  follows.  First,  the  BNP  calibration  data  taken 
from  a  pure  BNP  solution  at  various  oxygen  levels  (described 
above)  were  preprocessed  by  normalizing  them  to  their  mean 
intensity.  Next,  they  were  spectrally  decomposed  using  a  two- 
component  non-negative  PARAFAC  algorithm.19  This  yielded 
the  component  spectra  of  the  fluorescence  and  phosphorescence 
emission,  which  can  be  seen  in  Fig.  1(a).  The  measured  spectra 
of  BNP  emission  can  then  be  expressed  as  a  linear  combina¬ 
tion  of  each  of  the  two  components  using  a  non-negative  least- 
squares  fitting  algorithm.20  This  provides  the  F  and  P  component 
emission  intensities,  from  which  the  F/P  ratio  was  calculated. 
Finally,  this  was  calibrated  to  absolute  concentration  using  a  sec¬ 
ond  order  Stern  Volmer  model,  F/P  =  a  +  b[ 02]  +  c[ 02]2.  Fig¬ 
ure  1(b)  shows  the  results  of  this  calibration,  with  the  true  oxygen 
concentration  on  the  x-axis  and  the  extracted  oxygen  concen¬ 
tration  on  the  y-axis,  shown  as  a  boxplot.  It  was  found  that  the 
responsiveness  of  the  BNPs  plateaus  at  ~15%  oxygen,  with  bet¬ 
ter  accuracy  at  lower  oxygen  levels.  The  root- mean-square  error 
on  an  independent  testing  data  set  was  found  to  be  2.1  over  the 
entire  range,  or  1.5  when  including  only  oxygen  levels  of  15%. 

1 .4  Window  Chamber  Model 

In  vivo  imaging  was  performed  using  a  dorsal-skin-fold  window- 
chamber  model  in  a  nude  mouse  (source,  NCI  Frederick).21 
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Wavalength  (nm) 


True  Oxygenation  (%) 

Fig.  1  Component  spectra  (a)  are  shown  for  fluorescence  and  phos¬ 
phorescence  emission.  In  addition,  (b)  the  accuracy  of  extraction  of 
the  percent  of  oxygen  is  shown  for  the  calibration  data,  where  it  can 
be  seen  that  there  is  tight  correlation  at  lower  values,  with  saturation 
occurring  at  ~1 3%  oxygen. 


Briefly,  this  consists  of  a  titanium  metal  frame  that  is  surgically 
implanted  onto  a  dorsal  skin  fold.  The  frame  holds  the  skin 
taught,  and  one  face  of  the  skin  is  surgically  excised,  exposing 
the  underlying  skin  tissue,  an  example  of  which  is  shown  in 
Fig.  2.  Approximately  10,000  4T1  mammary  carcinoma  tu¬ 
mor  cells  were  injected  in  20  /xL  of  Debucco’s  modified  eagle 
medium  (DMEM)  via  a  30-gauge  syringe.  These  cells  expressed 
constitutive  red  fluorescent  protein  (RFP),  with  HIF-1  inducible 
GFP  expression.18  Finally,  a  coverglass  was  placed  over  the  tu¬ 
mor  and  secured  with  a  retaining  ring.  The  data  presented  were 
acquired  from  a  single  mouse.  All  animal  work  was  approved 
by  Duke  University’s  Institutional  Animal  Care  and  Use  Com¬ 
mittee. 


1 .5  In  Vivo  Imaging 

The  mouse  was  anesthetized  using  1.5%  isoflurane  mixed  with 
medical  air.  The  animal’s  temperature  was  maintained  using 
a  heated  blanket.  Breathing  was  monitored  during  imaging  to 
ensure  appropriate  depth  of  anesthesia.  First,  a  set  of  background 
images  were  acquired  to  characterize  tissue  autofluorescence. 
Next,  ~100  /xL  of  the  BNP  solution  was  injected  into  the  space 
between  the  window  and  tumor  tissue  to  bathe  the  tissue,  with 
excess  fluid  draining  out  the  sides  of  the  window,  leaving  a  thin 
layer  of  BNP  solution.  Fluorescence  images  were  acquired  every 

7.5  min  for  1  h,  with  transmission  images  acquired  continuously 
in-between  fluorescence  acquisitions. 


Fig.  2  The  window-chamber  model  used  in  this  study  consists  of  a 
titanium  metal  frame,  which  holds  the  dorsal  skin  fold  in  place.  The 
front  skin  flap  is  excised,  and  tumor  cells  are  injected  into  the  back 
plane  of  the  opposing  side  skin  and  covered  by  a  round  coverglass. 
The  inset  shows  the  mouse  with  window  implanted,  with  a  larger  view 
of  the  window  itself.  This  is  a  different  animal  than  that  shown  in  the 
other  figures. 

1.6  Data  Processing 

Transmission  images  were  processed  using  a  modified  form  of 
existing  techniques.22  A  modified  form  of  Beer’s  law  was  used, 
which  approximates  attenuation  through  the  tissue  as  a  linear 
combination  of  absorption  and  scattering  terms,  namely 

lo8io  (y)  =  b0  +  b\HeK  +  Y^s‘ci’  (!) 

'  '  i 

where  Iq  is  the  transmission  in  the  absence  of  tissue  attenuation, 
which  was  measured  using  neutral  density  filters;  b o  is  a  constant 
term  that  accounts  for  variability  in  source  intensity,  specular  re¬ 
flection,  or  other  heterogeneities  in  overall  intensity;  and  /xe &  is 
the  effective  attenuation  coefficient23  derived  from  diffusion  the¬ 
ory,  /xeff  =  y3/xfl(/xa  +  /x'),  calculated  for  representative  skin 
optical  properties  at  each  wavelength,24  and  accounts  for  attenu¬ 
ation  due  to  nonhemoglobin  absorption  and  scattering.  This  was 
empirically  found  to  provide  a  good  fit  to  the  nonvascular  com¬ 
ponent  of  attenuation  in  tissue  (for  example,  in  the  nonvascular 
regions  of  the  images).  st  is  the  wavelength-dependent  molar 
extinction  coefficient  of  each  absorber  in  the  model,  and  Q  is 
the  product  of  its  concentration  and  the  effective  path  length.  In 
this  case,  two  absorbers,  oxygenated  (Hb02)  and  deoxygenated 
(Hb)  hemoglobin,  were  used.  The  use  of  this  approximation 
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Fig.  3  (a)  Bright  field,  (b)  oxygen  tension  in  the  percent  of  oxygen,  (c)  hemoglobin  oxygen  saturation,  (d)  GFP-HIF-1  fluorescence,  and  (e)  consti- 
tutively  expressed  RFP  fluorescence  are  shown.  The  hemoglobin  saturation  (c)  is  shown  with  percent  oxygen  saturation  on  a  color  scale,  whose 
brightness  is  modulated  by  the  total  hemoglobin  content  (thus,  well-vascularized  regions  appear  brighter).  Also  in  (f),  a  combined  plot  shows  the 
hemoglobin  oxygen  saturation  plotted  as  in  (c),  with  a  3-D  projection  of  the  tissue  oxygenation  shown  on  the  vertical  axis.  The  tumor  is  seen 
predominantly  within  the  black  box.  These  images  were  acquired  toward  the  end  of  the  imaging  sequence  (t  =  55  min).  Images  are  approximately 
5  x  3.7  mm,  with  a  1-mm  scale  bar  seen  in  (a). 


enables  rapid  calculation  of  physiologic  parameters,  namely, 
hemoglobin  content  and  saturation  on  a  pixel-by-pixel  basis. 
This  equation  is  fit  using  a  non-negative  least-squares  optimiza¬ 
tion  algorithm  over  the  measured  wavelength  range.20  From 
this,  the  total  hemoglobin  content  times  the  path  length,  Hbtot  = 
CHbo2  +  CHb  and  the  hemoglobin  saturation,  Hbsat  =  CHbo2/ 
Hbtot  are  calculated. 

2  Results  and  Discussion 

Figure  3  provides  an  example  of  optical  data  demonstrating 
quantitative  imaging  of  tumor  hypoxia  and  vascular  function. 
Figure  3(a)  shows  a  transmission  image  that  allows  for  visual¬ 
ization  of  the  tumor  and  its  vasculature.  Figure  3(b)  depicts  the 
oxygen  tension  as  percent  O2  extracted  from  the  BNP  spectrum. 
These  images  reveal  a  region  in  the  tumor  with  relatively  high 
oxygen  tension,  while  relatively  low  oxygen  tension  appears  in 
the  adjacent  regions.  This  poorly  oxygenated  zone  appears  to  be 
poorly  vascularized  and  downstream  of  the  tumor  blood  supply. 
Figure  3(c)  shows  the  hemoglobin  saturation  on  a  color  scale 
from  blue  (deoxygenated)  to  red  (oxygenated).  The  color  en¬ 
coding  at  each  pixel  is  then  multiplied  by  the  total  hemoglobin 
content  extracted  from  the  same  pixel,  such  that  large  vessels 
appear  brighter,  while  nonvascularized  pixels  appear  black,  due 


to  low  hemoglobin  content.  The  tumor  appears  well  vascular¬ 
ized  with  heterogeneous  oxygenation.  Figure  3(d)  shows  the 
GFP  intensity,  normalized  to  the  peak  value,  which  shows  dis¬ 
tinct  pockets  of  intense  GFP  fluorescence  within  the  central 
part  of  the  tumor.  Figure  3(e)  shows  the  constitutively  active 
RFP,  showing  the  extent  of  the  tumor.  Finally,  Fig.  3(f)  shows  a 
three-dimensional  (3-D)  visualization  of  the  relationship  among 
hemoglobin  oxygen  saturation,  hemoglobin  content,  and  tissue 
oxygen  tension.  The  hemoglobin  oxygen  saturation  color  map  is 
overlaid  on  a  3-D  projection  of  the  tissue  oxygen  tension  plotted 
on  the  surface  map,  with  the  oxygen  tension  appearing  on  the 
vertical  axis.  Thus,  oxygen  gradients  can  be  visualized  as  slopes 
on  the  plot,  with  a  peak  appearing  on  one  side  of  the  tumor  close 
to  the  large  artery  on  the  lower  left  side,  with  a  steep  fall  off 
on  the  trailing  edge.  Video  1  shows  this  rendering  as  it  appears 
from  a  variety  of  angles. 

A  key  advantage  of  both  BNP  emission  and  hemoglobin 
absorption-based  hyperspectral  imaging  techniques  is  their  abil¬ 
ity  to  monitor  dynamic  changes  over  time.  To  this  end,  the  pa¬ 
rameters  described  for  Fig.  3  were  imaged  semicontinuously  for 
the  course  of  1  h.  This  enabled  characterization  the  temporal 
dynamics  of  each  parameter  and  allowed  for  their  interrelation¬ 
ships  to  be  studied.  Figure  4(a)  shows  the  Spearman  rank  cor¬ 
relation  between  the  hemoglobin  saturation  and  tissue  oxygen 
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Video  1  The  3-D  projection  of  tumor  oxygen  tension  on  the  z-axis, 
along  with  hemoglobin  saturation  and  total  hemoglobin  on  the  color 
scale  identical  to  that  shown  in  Fig.  3  (blue  =  deoxygenated  ->  red  = 
oxygenated,  with  the  intensity  modulated  by  the  total  hemoglobin  con¬ 
tent  making  vessels  appear  brighter).  This  enables  visualization  of  gradi¬ 
ents  in  tissue  oxygen  tension  in  relation  to  vessel  morphology  and  oxy¬ 
genation.  It  can  be  seen  that  the  tumor  oxygen  tension  is  highly  hetero¬ 
geneous  and  has  a  steep  gradient  dividing  a  relatively  well -oxygenated 
region  from  a  relatively  hypoxic  region.  (QuickTime,  2.6  MB). 

[URL:  http://dx.d0i.0rg/l 0.1 1 1  7/1 .3523363.1  ] 


Video  2  Shows  the  dynamics  of  hemoglobin  saturation  and  tumor  oxy¬ 
genation  over  an  1  -h  time  period  with  oxygen  tension  shown  as  a  3-D 
projection  along  the  z-axis,  and  hemoglobin  saturation  shown  on  a 
blue  red  color  scale  identical  to  that  in  Fig.  3.  Tissue  oxygenation 
can  be  seen  to  be  a  dynamic  process  with  constant  variability  through¬ 
out  the  time  period.  Hemoglobin  saturation  tracks  well  with  changes  in 
tissue  oxygenation,  with  reddening  vessels  coinciding  with  increased 
oxygen  tension.  (QuickTime,  3.4  MB). 

[URL:  http://dx.d0i.0rg/l  0.1 1 1  7/1 .3523363.2] 


tension,  pixel  by  pixel.  Solid  white  regions  indicate  areas  where 
the  hemoglobin  content  was  not  sufficient  to  reconstruct  the 
hemoglobin  saturation.  It  can  be  seen  that  there  is  a  tight  corre¬ 
lation  of  hemoglobin  oxygen  saturation  with  tumor  oxygenation 
(median  r  =  0.77).  This  correlation  is  relatively  low;  however, 
in  the  central  tumor  region,  where  instable  flow  and  shunting  is 
expected  to  occur.  Figure  4(b)  shows  the  pixel-by  pixel  standard 
deviation  of  tissue  oxygenation  over  time,  which  indicates  the 
magnitude  of  oxygen  fluctuations  over  time.  It  can  be  seen  that 
there  are  two  distinct  regions  within  the  tumor,  with  the  lower 
left  quadrant  showing  distinctly  greater  fluctuations.  This  region 
is  sandwiched  between  the  large  artery  and  vein  in  the  lower  left 
quadrant.  The  region  to  the  right  of  the  tumor  shows  consistently 
lower  oxygen  levels  throughout  the  imaging  period,  and  thus  has 
relatively  low  standard  deviation.  Video  2  shows  the  dynamics 
of  tissue  oxygen  tension  and  hemoglobin  saturation  over  time, 
which  illustrates  these  trends  visually.  Qualitatively,  it  can  be 
seen  that  there  are  often  rapid  changes  in  hemoglobin  saturation 


(color  scale),  that  coincide  with  more  gradual  changes  in  tissue 
oxygenation  (z-axis  height). 

Finally,  insofar  as  it  is  widely  hypothesized  that  different 
microenvironmental  phenotypes  respond  differently  to  various 
forms  of  treatment,  one  of  the  primary  aims  of  this  work  was 
to  evaluate  the  ability  of  these  imaging  techniques  to  segment 
tissues  into  different  physiologic  phenotypes  based  on  the  tem¬ 
poral  dynamics  of  tissue  oxygenation.  To  this  end,  the  temporal 
oxygenation  data  were  first  decomposed  using  wavelet  analysis. 
This  is  similar  to  Fourier  analysis,  which  has  previously  been 
used  to  characterize  temporal  dynamics,25  but  has  the  advantage 
of  being  localized  with  respect  to  time  as  well  as  frequency,26 
and  thus  would  better  describe  potentially  abrupt  or  irregular 
changes  in  flow  patterns.  A  decomposition  was  performed  using 
a  second-order  Daubechies  wavelet 26  The  discrete  wavelet  de¬ 
composition  is  applied  at  multiple  levels,  or  scales.  At  the  first 
level,  the  original  signal  is  decomposed  into  an  approximation 
and  detail  signal,  where  the  detail  signal  is  the  convolution  of 
the  original  signal  and  the  wavelet  function  at  the  smallest  scale 
or  highest  frequency.  This  essentially  splits  the  signal  into  a 


Fig.  4  (a)  The  Spearman  rank  correlation  of  hemoglobin  oxygen  saturation  and  tissue  oxygen  tension,  which  indicates  that  the  central  portion  of 
the  tumor  has  relatively  low  correlation  between  these  parameters,  indicating  inefficiency  of  oxygen  delivery  to  the  tumor  and  possible  shunting  of 
flow,  (b)  The  standard  deviation  of  the  tissue  oxygen  tension  over  the  1  -h  time  period  over  which  measurements  were  made,  which  indicates  regions 
of  the  tumor  that  show  the  largest  magnitude  of  fluctuations  over  time.  A  1  -mm  scale  bar  is  seen  in  (a). 
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Fig.  5  Results  of  the  clustering  analysis  shows  (a)  the  spatial  cluster  assignment  and  (b)  the  median  percent  O2  within  each  cluster.  Also  shown  is 
the  (c)  HIF-1 -induced  GFP  expression  for  the  same  tissue  region.  There  is  clear  spatial  segregation  of  the  tumor  and  normal  tissues,  with  the  tumor 
being  contained  predominantly  within  clusters  2,  3,  and  4.  Note  that  this  analysis  was  performed  over  a  smaller  region  of  interest  indicated  by  the 
black  box  in  Fig.  3(a).  A  1-mm  scale  bar  is  seen  in  (a). 


high-frequency  component  (detail)  and  low-frequency  compo¬ 
nent  (approximation),  the  sum  of  which  yields  the  original  sig¬ 
nal.  This  process  can  be  repeated  with  the  scale  of  the  wavelet 
function  increasing  by  powers  of  2  at  each  step,  thus  extracting 
lower  and  lower  frequency  components.  In  this  case,  a  two-level 
decomposition  was  performed,  which  divides  the  temporal  re¬ 
sponse  into  short-  and  long-period  oscillations  (center  cycling 
periods  of  approximately  22.5  and  45  min).  These  coefficients 
were  fed  into  a  fuzzy  c-means  clustering  algorithm27  to  identify 
which  regions  within  and  surrounding  the  tumor  behaved  most 
similarly  over  time.  The  fuzzy  clustering  algorithm  classifies 
each  tissue  region  according  to  the  relative  associations  with 
each  cluster  centroid,  allowing  for  the  fact  there  is  not  a  hard 
discrimination  boundary  between  clusters;  thus,  cluster  associ¬ 
ations  are  assigned  weights  from  0  to  1.  This  was  performed 
over  a  smaller  region  of  interest  containing  only  the  tumor  and 
immediately  adjacent  normal  tissue  as  shown  by  the  dark  box 
in  Fig.  3(a)  to  simplify  this  analysis.  For  the  purpose  of  display, 
the  maximum  cluster  association  was  selected  for  each  pixel. 
Figure  5(a)  shows  a  map  of  the  maximum  cluster  association 
for  each  pixel,  and  Fig.  5(b)  shows  the  median  percent  O2  over 
time  within  each  of  these  clusters.  The  tumor  is  contained  pre¬ 
dominantly  within  clusters  2,  3,  and  4,  with  cluster  3  showing 
an  earlier  increase  in  oxygenation  [see  Fig.  5(b)]  and  also  hav¬ 
ing  a  higher  median  value  and  larger  variance  over  time  [see 
Fig.  3(b)].  Interestingly,  the  bimodal  distribution  of  GFP  ex¬ 
pression  is  also  reflected  in  the  distribution  of  clusters  3  and  4, 
which  may  reflect  the  influence  of  oxygen  dynamics  on  HIF- 
1  activity.  Specifically,  cluster  3  has  the  largest  fluctuations  in 
oxygen  tension  [as  seen  in  Figs  4(b)  and  5(b),  whereas  cluster 
4  has  the  lowest  absolute  oxygen  tension  [as  seen  in  Fig.  5(b)], 
which  have  been  shown  to  be  two  significant  drivers  of  HIF-1 
activity.4  A  pixel-by-pixel  correlation  of  cluster  3  with  GFP  fluo¬ 
rescence  intensity  yields  a  Pearson  correlation  coefficient  of  0.46 


(p  <  lx  10-15),  indicating  the  significant  correlation  between 
cycling  hypoxia  dynamics  and  HIF-1  activity.  The  remaining 
clusters  were  not  correlated  with  GFP  fluorescence  (correlation 
coefficients  ranging  between  —0.14  and  —0.09),  which  may 
be  due  in  part  to  the  fact  that  they  all  include  significant  normal 
tissue  regions. 

One  of  the  limitations  of  this  study  is  the  Beer’s  law  assump¬ 
tion  for  tissue  attenuation.  This  is  a  common  practice  in  inter¬ 
preting  wide-field  hyperspectral  imaging  data  sets,  but  it  does 
not  account  for  changes  in  path  length  with  respect  to  wave¬ 
length  and  optical  properties,  or  for  the  heterogeneity  present 
in  tissue.  We  are  currently  developing  a  Monte  Carlo-based 
approach  for  interpreting  such  data  and  evaluating  the  errors  in¬ 
duced  through  use  of  simplified  modeling  techniques.  Another 
consideration  in  interpreting  these  data  is  the  route  of  admin¬ 
istration.  Because  the  nanoparticles  were  applied  topically,  the 
sensing  volume  is  limited  to  the  superficial  tissue  space.  It  is 
unclear,  based  on  this  study,  whether  the  nanoparticles  are  taken 
up  by  the  cells  or  to  what  extent  they  penetrate  into  the  tissue, 
but  this  could  influence  the  measurement  over  time  as  the  sens¬ 
ing  volume  changes.  Because  epi-illumination  was  used,  along 
with  a  relatively  short  imaging  duration  (1  h),  it  is  likely  that  the 
measurement  is  highly  superficial  throughout  the  measurement 
time.  However,  these  fluctuations  are  more  likely  driven  by  true 
changes  in  tissue  oxygenation  than  a  change  in  penetration  depth 
over  time,  given  the  high  correlation  coefficient  seen  between 
the  independent  measure  of  hemoglobin  oxygen  saturation.  The 
time  scale  for  these  fluctuations  in  oxygenation  is  consistent 
with  those  measured  previously  using  other  methods.4,25  In  the 
future,  confocal  imaging  or  histology  could  be  employed  to  gain 
a  better  understanding  of  the  penetration  depth  of  these  nanopar¬ 
ticles  over  time,  after  topical  application.  In  addition,  for  many 
clinical  applications,  a  topical  application  may  not  be  practical. 
For  this  reason,  the  use  of  intravenous  administration  has  also 
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been  considered  and  we  have  developed  a  PEGylated  nanopar¬ 
ticle  formulation  that  was  found  to  accumulate  near  the  tumor 
vasculature  using  confocal  imaging.28 

In  conclusion,  we  have  demonstrated  and  validated  tech¬ 
niques  by  which  vascular  and  tissue  oxygenation  can  be  imaged 
at  high  spatial  and  temporal  resolution,  enabling  a  detailed  un¬ 
derstanding  of  the  tumor  microenvironment  and  its  influence  on 
tumor  gene  expression.  All  these  techniques  are  implemented 
using  a  simple  modification  to  a  standard  fluorescent  micro¬ 
scope,  enabling  a  wide  range  of  potential  applications.  There  is 
no  other  technology  that  facilitates  simultaneous  measurement 
of  vascular  oxygen  supply  and  tissue  oxygen  tension  using  a  sin¬ 
gle  instrument,  and  this  work  enables  a  significant  advancement 
in  our  ability  to  characterize  tissue  hypoxia  and  its  dynamic  in¬ 
teractions  in  vivo.  Hypoxia  is  a  critical  factor  in  a  wide  range 
of  disease  processes,  including  cancer,  epilepsy,  wound  healing, 
and  more.  The  capability  of  the  techniques  developed  here  to 
quickly,  quantitatively,  and  nondestructively  assess  the  dynam¬ 
ics  of  vascular  oxygen  supply  and  tissue  oxygen  tension  could 
open  up  new  possibilities  for  understanding  of  the  basic  mech¬ 
anisms  by  which  disease  processes  are  influenced  by  hypoxia 
and  correlating  these  to  therapeutic  response.  In  the  future,  we 
hope  to  translate  these  techniques  to  a  clinical  setting  as  a  means 
of  assessing  tumor  hypoxia  and  therapeutic  response  to  evaluate 
their  potential  in  guiding  and  monitoring  therapy. 
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